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ABSTRACT 

Modern  naval  warships  rely  on  vast  arrays  of  sensor  networks  to  evaluate  the  performance  of 
mission  critical  systems.  Although  these  sensor  networks  enable  increased  levels  of  automation, 
they  are  costly  to  install  and  to  maintain.  The  power  distribution  network  offers  an  alternative 
solution  for  tracking  the  performance  of  mission  critical  systems.  Research  conducted  at 
Massachusetts  Institute  of  Technology’s  Laboratory  for  Electromagnetic  and  Electronic  Systems 
(LEES)  has  proven  that  the  power  distribution  network  contains  vital  information  that  can 
provide  performance  monitoring  and  automatic  diagnostic  functions. 

This  thesis  will  address  the  issue  of  sensor-count  reduction  through  the  application  of  Non- 
Intrusive  Load  Monitoring  (NILM)  technology.  Theoretical  studies  and  field  experiments  will 
be  presented  in  order  to  demonstrate  the  NILM’s  ability  to  correlate  load  activity  with  power 
measured  from  an  aggregate  level  in  the  distribution  system.  Additionally,  a  critical  evaluation 
is  conducted  on  the  current  NILM  configuration’s  ability  to  perform  automated  classification. 
Findings  will  be  supported  using  data  collected  from  NILMs  monitoring  power  flow  on  board 
the  U.S.  Coast  Guard  Cutter  ESCANABA  (WMEC-907). 
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Chapter  1 


Introduction 
1. 1  Background 

The  Non-Intrusive  Load  Monitor  (NILM)  is  a  device  used  to  measure  raw  current  and  voltage 
from  a  single  point  in  an  electrical  distribution  system  (Figure  1-1).  These  measurements  are 
converted  (preprocessed)  into  real  and  reactive  power  (Figure  1-2)  and  their  associated 
harmonics.  When  applied  to  a  given  electro-mechanical  system,  the  power  measurements 
provide  vital  system  information  that  can  be  used  to  determine  the  overall  “health”  of  not  only 
the  electrical  components  (i.e.  motors),  but  also  the  mechanical  components  (i.e.  pumps,  filters, 
valves,  etc.)  of  a  given  system. 


Time  <s} 

Figure  1-1:  Raw  AC  voltage  and  current  measurements  recorded  during  a  motor  start-up.  (1) 


Time  (s) 

Figure  1-2:  Preprocessed  spectral  power  envelopes  of  motor  transient  in  Figure  1-1.  The  high  frequency  content 
associated  with  the  raw  data  has  been  removed.  Preprocessor  information  is  typically  used  for  NILM  applications  that 

involve  real-time  detection  and  classification  algorithms.  (1) 
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The  NILM’s  effectiveness  in  shipboard  power-based  evaluation  and  monitoring  has  been 
demonstrated  on  several  occasions.  Early  research  focused  on  collecting  power  data  from 
individual  loads  (auxiliary  seawater  system,  sewage  system,  reverse  osmosis),  then  analyzing  the 
archived  data  to  determine  correlations  between  system  operation  and  power  signatures.  (2)  (3) 
(4)  More  recently,  research  efforts  have  focused  on  applying  the  knowledge  gained  from 
individual  load  monitoring  to  modern  state  and  parameter  estimation  algorithms  to  provide 
system  “health”  information  for  user  interfaces  and  diagnostic  tools.  (5)  (6)  (7) 

Figure  1-3  shows  a  NILM  layout,  where  multiple  loads  are  monitored  via  a  single  NILM.  As  the 
NILM  moves  further  “upstream”  to  monitor  multiple  loads  from  a  single  point,  it  becomes  less 
“intrusive”.  However,  there  are  practical  limitations  to  how  far  “upstream”  the  NILM  can 
effectively  perform  such  tasks  as  event  detection  and  system  diagnostics.  From  a  designer’s 
perspective,  there  are  necessary  tradeoffs  that  impact  the  NILM’s  performance  in  real  world 
applications.  This  thesis  will  address  the  tradeoffs  that  affect  the  likely  performance  of  the 
NILM  in  shipboard  power  systems.  Additionally,  current  event  detection  and  classification 
methods  will  be  scrutinized  using  both  local  and  upstream  data  collected  during  recent  field 
studies  on  the  USCGC  ESCANABA  (WMEC-907)  and  the  U.S.  Navy’s  Land  Based 
Engineering  Site  (LBES)  in  Philadelphia,  Pennsylvania. 


Figure  1-3:  Typical  system  layout  for  a  NILM  used  to  monitor  multiple  loads.  (1) 
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1.2  Motivation  for  Research 


The  U.S.  Department  of  the  Navy  spends  approximately  $4  billion  a  year  to  maintain  its  existing 
fleet  (8).  In  an  effort  to  lower  the  costs  associated  with  maintaining  an  aging  fleet,  the  navy 
enterprise  has  incorporated  policy  aimed  at  reducing  lifecycle  maintenance  costs.  In  May  of 
1998,  the  Chief  of  Naval  Operations  instituted  a  condition-based  maintenance  (CBM)  policy  (9), 
which  focused  on  eliminating  waste  associated  with  traditional  navy  maintenance  practices  that 
relied  on  preventative  maintenance  strategies. 

During  the  same  time  period  as  CBM  implementation,  the  Smart  Ship  concept  was  also  gaining 
momentum  (10).  The  focus  of  the  Smart  Ship  project  was  to  identify  key  enablers  for  reducing 
manning  and  workload.  Initially,  61  enabling  technologies  were  identified  and  installed  on  the 
navy’s  prototype  Smart  Ship,  the  USS  YORKTOWN  (CG  48).  These  enablers  were  divided  into 
seven  core  technologies,  one  of  which  was  the  Integrated  Condition  Assessment  System  (ICAS) 
(11). 

ICAS  is  the  U.S.  Navy’s  CBM  program  of  record;  it  is  currently  installed  on  95  surface  ships 
(10).  The  ICAS  provides  the  capability  to  collect,  monitor,  and  trend  machinery  data  via  various 
types  of  sensors  located  within  the  systems  of  interest.  This  information  can  then  be  used  locally 
by  ship’s  force  to  diagnose  equipment  or  sent  to  intermediate  maintenance  facilities  for  diagnosis 
and  trend  analysis.  However,  in  order  for  ICAS  to  realize  its  full  potential,  it  must  continue  to 
grow  and  increase  its  capabilities  to  meet  the  needs  of  the  future  navy. 

One  area  for  potential  improvement  in  ICAS  technology  is  sensor  count  reduction.  ICAS  relies 
on  an  extensive  array  of  sensors  to  provide  the  vital  information  regarding  component  status.  As 
the  sensor  count  grows,  so  does  the  complexity  and  costs  associated  with  maintaining  the  very 
system  designed  to  reduce  maintenance  costs.  The  USS  San  Antonio  (LPD-17)  has  7023  signals 
associated  with  the  Machinery  Control  System  (MCS),  while  the  amphibious  assault  ships  have 
over  12,000  signals.  The  Navy’s  Future  Surface  Combatant  (DDG  1000)  is  projected  to  have  up 
to  33,000  control  signals.  (12) 
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For  a  little  more  than  a  decade,  researchers  at  the  Massachusetts  Institute  of  Technology  (MIT) 
have  analyzed  ways  to  monitor  shipboard  systems  using  the  Non-Intrusive  Load  Monitor 
(NILM).  The  NILM  is  a  system  that  can  determine  the  operating  schedule  of  electrical  loads  in  a 
target  system  using  centralized  measurements  (13).  By  monitoring  the  power  signatures 
associated  with  individual  loads  from  an  electrical  node,  the  NILM  has  the  potential  to 
drastically  reduce  sensor  count  currently  associated  with  U.S.  Navy  vessels.  This  sensor  count 
reduction  combined  with  the  existing  ICAS  architecture  could  prove  tremendously  beneficial  to 
the  Navy’s  condition  based  maintenance  program. 

1.3  Objective  and  Outline 

Previous  shipboard  research  involving  the  NILM  relied  on  dedicated  component  monitoring. 
This  approach  allowed  researchers  to  effectively  demonstrate  several  basic  diagnostic  indicators 
for  mechanical  loads  based  on  electrical  power  signatures.  Current  research  has  focused  on 
using  the  NILM  to  conjointly  monitor  several  loads  from  an  aggregate  point  in  the  electrical 
power  system.  The  objective  of  this  research  is  to  explore  how  well  the  NILM  can  determine 
load  composition  within  shipboard  aggregate  power  systems  using  spectral  envelope  estimates. 

The  thesis  consists  of  five  chapters  including  the  introduction  and  conclusions.  Chapter  2  takes 
an  in-depth  look  at  three  key  challenges  that  must  be  addressed  when  designing  aggregate  power 
monitoring  systems.  Methods  for  aggregate  event  detection  that  are  identified  in  Chapter  2  are 
applied  in  Chapter  3,  which  uses  data  gathered  from  the  Navy’s  Land  Based  Engineering  Site 
(LBES)  in  Philadelphia,  Pennsylvania.  Chapter  4  provides  a  brief  description  of  the  most  current 
NILM  configuration,  which  uses  a  software  package,  Ginzu,  to  perform  automated  classification 
and  diagnostics.  Lessons  learned  from  a  recent  deployment  of  the  USCGC  ESCANABA 
(WMEC-907)  with  this  particular  NILM  installed  are  also  explored. 
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Chapter  2 


Aggregate  Power  System  Monitoring 


The  power  distribution  network  can  be  exploited  not  only  as  a  means  for  power  distribution,  but 
also  as  a  means  for  providing  useful  diagnostic  system  monitoring.  This  chapter  focuses  on 
recent  research  conducted  to  demonstrate  how  well  the  NILM  performs  in  an  aggregate 
environment,  where  the  objective  is  to  minimize  the  number  of  NILM  monitors  required  to 
capture  useful  system  information  for  diagnostic  monitoring.  Figures  and  equations  provided 
throughout  this  chapter  were  taken  from  (1),  which  was  co-authored  by  the  author  of  this  thesis. 

Ideally,  a  single  NILM  would  be  able  to  perform  diagnostic  monitoring  of  an  entire  ship; 
however,  there  are  practical  limitations  that  prevent  this  level  of  “non-intrusiveness”.  Although 
not  all  inclusive,  the  following  three  areas  heavily  influence  the  NILM’s  ability  to  perform 
aggregate  diagnostic  monitoring  and  are  the  focal  point  of  this  chapter: 

•  First,  the  NILM  samples  a  continuous  signal  (i.e.  current),  which  results  in  quantization. 
This  quantization  is  the  result  of  the  analog- to-digital  conversion  (ADC),  where  the 
analog  signal  must  be  converted  into  finite  digital  data  output.  For  a  given  NILM 
hardware  configuration,  the  power  range  it  observes  is  bounded  by  the  sensor  and  ADC 
card  selection.  As  the  NILM  is  moved  further  “upstream”  in  the  power  distribution 
system,  the  corresponding  power  range  increases.  The  larger  power  range  results  in 
different  hardware  selection,  which  ultimately  influences  the  NILM’s  ability  to 
accurately  detect  smaller  transients. 

•  Second,  the  NILM  is  susceptible  to  external  signals  that  distort  underlying  information  in 
acquired  signals.  This  includes  random  measurement  noise  and  other  distortions  such  as 
possible  electromagnetic  noise. 

•  Third,  the  NILM  must  be  capable  of  disaggregating  loads  of  interest.  For  example,  there 
may  be  two  induction  motors  of  similar  size  associated  with  separate  systems.  Detection 
techniques  that  rely  on  power  level  changes  alone  may  not  be  capable  of  distinguishing 
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the  two  motors.  Additionally,  random  power  changes  due  to  such  loads  as  cycling  lube 
oil  heaters  may  mask  transients  of  interest  for  a  given  detection  algorithm. 

2.1  NILM  System 

Figure  2-1  illustrates  how  continuous  signals  (current  and  voltage)  are  (1)  sensed,  (2)  sampled, 
and  then  (3)  processed  to  perform  various  load  monitoring  functions.  A  typical  NILM 
configuration  consists  of  an  enclosed  sensor  box,  external  power  connections  to  include  current 
transducers  and  voltage  sensing  lines,  and  a  personal  laptop  computer  for  data  processing  and 
archiving. 


Sensing 


Processing 


Figure  2-1:  NILM  data  acquisition  and  processing  diagram. 


Acquired  data  in  the  form  of  current  and  voltage  is  first  converted  into  digital  form  through  an 
analog-to-digital  converter.  The  quantized  signal  from  the  converter  is  passed  to  a  personal 
computer  (PC)  via  Ethernet  for  preprocessing.  Preprocessed  data  extracts  useful  information 
from  the  raw  data  signal  for  detection,  classification,  and  diagnostic  functions.  Figure  2-2 
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provides  an  internal  view  of  a  NILM  sensor  box.  Figure  2-3  displays  the  physical  layout  of  an 
entire  NILM  install. 


Figure  2-2:  Integrated  NILM  sensor  box  courtesy  of  NEMOmetrics.  Sensor  box  includes  connections  for  power,  current 
signals,  as  well  as  an  Ethernet  connection  to  connect  to  an  external  PC. 


Figure  2-3:  NILM  layout  for  installation  on  USCGC  ESCANABA  (WMEC-907). 
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2.1.1  Data  Acquisition 


Data  acquisition  measures  the  current  and  voltage  that  supplies  power  to  a  particular  load  or 
group  of  loads  within  a  system.  Although  voltage  and  current  are  both  sensed,  this  discussion 
will  focus  on  current  alone.  A  majority  of  the  information  that  is  used  for  NILM  applications  is 
derived  from  the  current  signal.  A  current  transducer  such  as  the  one  shown  in  Figure  2-4  is 
placed  around  a  particular  phase  to  sense  the  incoming  primary  current,  Ip.  The  Ip  is  converted 
into  a  smaller  signal  by  a  conversion  ratio  KN.  The  output  of  the  CT  is  referred  to  as  the 
secondary  current,  Is .  TABLE  2-1  lists  commonly  used  current  transducers  for  NILM 
applications. 


Figure  2-4:  Interior  view  of  Power  Panel  2-82-3.  The  LA-205-S  current  transducer  (CT)  and  three  voltage  sensing  lines 

are  annotated. 
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TABLE  2-1:  TYPICAL  CURRENT  TRANSDUCERS  AND  CONVERSION  RATIOS 


Model 

Maximum  Primary  Current 

Kn 

LA-55 

±  50  d 

1/1000 

LA- 150 

±  150  A 

1/2000 

LA-205 

±200  A 

1/2000 

LA-305 

±  300  A 

1/2500 

The  secondary  output  Is  is  fed  to  a  burden  resistor  RM  to  ground.  Rmi s  proportional  to  Ip  and  can 
be  represented  as  follows: 


I'm  —  Is  *  Rm  —  Ip  *  Kn  *  Rm 


2-1 


For  an  LA-55  CT,  a  50  Q  burden  resistor,  and  an  Ip  of  ±  50  A, 

VM  =  ±50  A  *  *  50  n  2-2 


*  ±2.5  V 


2-3 


The  continuous  amplitudes  from  the  incoming  signal, LM,  are  applied  to  an  ADC  card.  The  ADC 
samples  the  signal  at  a  rate  of  8000  samples  per  second.  The  process  of  converting  the 
continuous  signal  into  digital  data  with  finite  precision  results  in  quantization.  The  number  of 
quantization  levels  L  can  be  derived  based  on  the  number  of  significant  binary  bits  B  available 
for  a  particular  ADC  card.  For  instance,  if  B  =3  bits,  there  are  2s  quantization  levels  or  8 
discrete  levels.  The  effect  of  quantization  on  an  incoming  continuous  signal  is  shown  in  Figure 
2-5.  As  the  bit  level  increases,  the  digital  output  closely  resembles  that  of  the  incoming 
continuous  waveform. 
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Figure  2-5:  Illustration  of  quantization  error  introduced  during  the  signal  processing  through  the  ADC.  The  continuous 
function,  sin(x),  is  compared  to  a  2  bit  and  4  bit  discrete  output  from  an  ADC. 

In  (1),  a  detailed  discussion  on  the  various  ADC  cards  typically  used  for  NILM  applications  is 
provided.  Currently,  the  LabJack  UE9  Ethernet-based  card  is  the  standard  for  NILM 
installations.  This  card  provides  design  flexibility  by  offering  Ethernet  connection  with  a 
downstream  PC.  Its  predecessor,  the  PCI- 17 10,  required  the  downstream  computer  to  be  placed 
near  the  sensor  box  to  accommodate  its  PCI  card.  The  LabJack  UE9  has  the  following 
specifications: 

•  Voltage  range  (Vr)  of  ±  5  V 

•  Counts  output  range  of  0  to  65 520 1 

•  12  bits  (B) 

One  way  to  describe  the  accuracy  of  a  given  ADC  card  is  by  measured  voltage  per  least 
significant  bit  (LSB),  where 

v  -  |VrI 

Vm.lsb  ~  2-4 


1  The  software  associated  with  the  LabJack  outputs  16-bit  values;  however,  there  are  only  12  significant  bits.  The 
lower  4  bits  have  zeros  for  placeholders.  This  results  in  the  count  output  range  listed. 
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Using  the  Lab  Jack  and  a  12  bit  ADC  converter, 


_  ((+5  V)  -  (-5  V)) 

VM,LSB  ~  212 

«  2.441  mV /bit 

This  represents  the  maximum  accuracy  for  a  particular  ADC  over  its  entire  input  voltage  range. 
Note  that  in  applications,  the  full  Vr  may  not  be  achieved;  this  value  is  dependent  upon  the 
burden  resistor  and  current  transducer  that  is  chosen.  Recall  from  2-3  that  the  measured  range 
was  only±  2.5  V,  which  results  in  a  loss  of  half  of  the  available  Vr.  This  loss  in  resolution  can 
be  mitigated  by  careful  selection  of  burden  resistors,  CTs,  and  ADC  cards  for  a  given  NILM. 


2-5 

2-6 


Another  useful  way  to  describe  the  accuracy  of  an  ADC  card  is  to  determine  the  primary  current 
per  LSB,  Ir,lsb ■  Using  the  components  in  the  example  above, 


Ip.LSB 


Vm.lsb 

kn  *  rm 


2.441  mV /bit 
Ip’LSB  ~  1/1000  *  5011 


2-7 


2-8 


«  48.82  mA/bit 


2-9 


which  yields  an  associated  primary  current  value  in  amps/bit.  Data  resolution  may  also  be 
described  by  a  comparison  of  a  target  load  current  I  to  IPiLSb ■  The  number  of  least  significant 
bits  N  corresponding  to  the  target  load  can  be  computed  as 

/ 

N  =  -  2-10 

Ip, LSB 

The  effective  number  of  bits  Be  is  found  by 

N  =  2Be  2-11 

=  l°g2(^0  2-12 
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BE  —  1°§2 


I 


2-13 


P.LSB 


For  instance,  if  I  =  ±  5  A  for  the  same  ADC,  burden  resistor,  and  CT  as  described  for  Equation 
2-9, 


Be  = 


l°g2 


(+5  A  —  (—5  A))\ 

48.82  mA/bit  ) 


2-14 


«  7.68  bits 


2-15 


Fully  expanded,  Equation  2-13  yields 


BE  ~  1°§2 


'I  *  KN  *  Rm  *  2 

v  K\ 
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2.1.2  Spectral  Envelope  Preprocessor 

The  NILM  utilizes  a  preprocessor  to  extract  spectral  power  envelopes  from  incoming  raw 
voltage  and  current  waveforms.  “Preprocessing  exposes  energy  consumption  and  harmonic 
signatures  that  are  relatively  uniquely  associated  with  the  physical  task  performed  by  a  load.”  (1) 
Similar  to  radio  transmissions  using  amplitude  modulation,  the  target  information  is  carried  on  a 
60  Hz  (for  typical  shipboard  distribution  systems)  waveform  with  varying  amplitude.  The  carrier 
waveform  is  removed,  revealing  the  power  spectrum  and  its  vital  information.  The  resulting 
envelope  is  used  to  identify  loads  and  perform  diagnostics. 

Reference  (14)  described  that  the  current  signal  x(t  )could  be  decomposed  using  Fourier  integrals 
into  the  in-phase  spectral  envelopes  au  and  the  quadrature  spectral  envelopes  b/(: 

2 

ajk(0  =  j 

bk(t)  =  y 


ft-T 

I  x(t)  sin(/t&)r)  dr 

ft-T 

I  x(t)  cos(/co)t)  dr 

>t 


2-17 


2-18 
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where  the  averaging  period  T  is  one  or  more  periods  of  the  line  frequency,  and  the  sinfAzurjis 
phase-locked  to  the  voltage  measurement  using  a  Kalman  filter.  (15)  The  integer  k  corresponds 
to  the  harmonic  index.  For  NILM  applications,  k  =  1  corresponds  to  real  and  reactive  power 
terms.  The  a*  term  provides  real  power  P  when  scaled  based  on  the  front-end  data  acquisition 
configuration.  Similarly,  reactive  power  Q  is  computed  using  the  bu  term  and  the  appropriate 
scale  factors.  Figure  1-2  illustrates  a  preprocessor  envelope  for  an  induction  motor. 

2.1.3  Load  Classification  and  Diagnostics 

To  realize  the  full  potential  of  the  NILM  for  shipboard  applications,  there  must  be  intelligent 
software  solutions  that  are  capable  of  distracting,  classifying,  and  diagnosing  transients. 
Methods  for  identifying  when  a  transient  occurs  often  rely  on  steady  state  power  levels  and 
transient  power  signatures.  Typically,  transients  are  noted  when  a  load  turns  on  and  off. 

Diagnostics  require  accurate  classification  of  transients  within  a  monitored  system.  Once 
properly  identified,  algorithms  can  be  used  to  look  for  “unusual”  behavior.  Take  the  CHT 
system  for  example.  An  exemplar  or  power  “fingerprint”  may  be  used  to  detect  when  a  vacuum 
pump  is  turned  on.  A  diagnostic  algorithm  could  then  look  for  unusual  behavior,  such  as 
excessive  run  times  or  improper  power  levels.  Information  may  then  be  displayed  to  the  operator 
to  indicate  improper  operations. 

Matching  the  transient  signatures  relies  on  finding  unique  aspects  of  a  load’s  power  usage. 
When  scaling  to  larger  NILM  applications,  the  ability  to  discern  these  unique  features  is 
diminished  as  the  number  of  loads  increases.  Smarter  methods  for  extracting  individual  load 
information  may  be  necessary.  This  is  addressed  in  further  detail  in  Section  2.5. 

2.2  Measurement  Accuracy 

In  order  to  successfully  perform  load  identification  and  diagnostics  of  the  aggregate  power 
signal,  the  NILM  must  maintain  a  certain  degree  of  accuracy  with  respect  to  data  measurements. 
Recall  from  previous  discussion,  that  the  designer  bounds  the  NILM’s  power  range  when 
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choosing  burden  resistors,  CTs,  and  ADC  cards.  As  the  NILM  moves  further  “upstream”,  load 
currents  increase  resulting  in  a  loss  in  resolution  for  smaller  loads  in  the  aggregate.  One  method 
for  evaluating  the  NILM’s  performance  to  disaggregate  loads  in  the  “upstream”  is  to  observe 
how  well  it  detects  and  classifies  target  loads. 

Two  NILMs  were  installed  on  the  USCGC  ESCANABA  on  February  14,  2009  to  measure  the 
effectiveness  of  the  NILM  to  disaggregate  load  information.  A  “downstream”  NILM  was 
installed  at  the  CHT  control  panel,  while  an  “upstream”  NILM  was  placed  at  a  node  in  the 
electrical  distribution  system  that  feeds  the  CHT  control  panel  and  other  downstream  loads.  The 
objective  was  to  describe  how  well  the  NILM  disaggregates  and  reports  the  operations  of 
individual  electrical  loads  using  aggregate  current  at  the  service  entry  to  a  power  panel  that 
supports  multiple  loads. 

The  upstream  node  for  the  installation  was  power  panel  P  2-282-3  (Figure  2-6),  which  is  located 
in  Auxiliary  Machinery  Room  #1,  centerline,  against  bulkhead  82.  It  receives  power  directly 
from  the  450  VAC,  60  hertz,  3cp,  800  amps  IS  main  distribution  switchboard.  The  panel’s 
isolation  breaker  has  a  current  trip  set  point  of  150  amps  and  is  located  on  the  IS  distribution 
board.  Panel  loads  range  from  the  CHT  control  panel  (largest  load  1 .7  HP  vacuum  pump)  to  the 
40  HP  fin  stabilizer  motor.  The  complete  list  of  loads  along  with  ratings  is: 

•  Wastewater  Discharge  Pump  #1  Motor  (5  HP) 

•  Fin  Stabilizer  #2  Motor  (40  HP) 

•  Reverse  Osmosis  Disconnect  Switch  (15  kW) 

•  HVAC  Hot  Water  Circulation  Pump  #2  Motor  (1.5  HP) 

•  Potable  Water  Service  Pump  #1  Motor  (7.5  HP) 

•  Potable  Water  Service  Pump  #2  Motor  (7.5  HP) 

•  CHT  Control  Panel  (Two  1.5  HP  Vacuum  Pump  Motors  and  Two  2  HP  Discharge  Pump 
Motors) 
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Figure  2-6:  Interior  view  of  P-2-282-3  (left)  and  close-up  view  of  bottom  portion  of  panel  (right). 


The  downstream  NILM  was  placed  at  the  local  control  panel  of  the  CHT  system  (Figure  2-7). 
The  control  panel  receives  power  from  the  upstream  power  panel,  P-2-282-3.  The  primary 
purpose  of  the  control  panel  is  to  provide  various  control  signals  to  the  CHT  system,  as  well  as 
provide  power  to  the  system’s  vacuum  and  discharge  pumps.2  Three  CTs  were  installed  in  order 
to  obtain  differently  scaled  raw  data  locally.  The  three  channels  were  connected  to  differently 
sized  burden  resistors  to  create  the  variation  in  data  scaling. 


2  Further  explanation  on  the  control  panel  can  be  found  in  Chapter  4. 
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Figure  2-7:  Interior  view  of  CHT  Control  Panel  (left)  showing  three  CTs  connected  to  a  single  phase  of  line  current. 
Interior  view  of  NILM  enclosure  box  (right)  shows  the  physical  location  of  CT  inputs  from  the  control  panel. 


A  single  run  of  CHT  vacuum  pump  transients  were  captured  locally  on  all  three  channels,  as  well 
as  “upstream”.  Although  three  different  local  configurations  were  installed,  the  following 
discussion  will  focus  only  on  the  downstream  configuration  with  the  highest  bit  resolution. 
TABLE  2-2  lists  the  NILM  front-end  components  that  provide  the  greatest  disparity  between 
“upstream”  and  “downstream”  scaling. 


TABLE  2-2:  FIELD  EXPERIMENT  NILM  COMPARISON 


Downstream 

Upstream 

Current  Transducer 

LA-55 

LA-205 

kn 

1/1000 

1/2000 

Rm 

155Q 

5Q 

Vr 

±sv 

±sv 

B 

12  bits 

12  bits 
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2.2.1  Raw  Data 


The  degree  of  closeness  with  which  the  “upstream”  NILM  can  measure  a  “downstream”  load  is 
diminished  due  to  selection  of  the  current  transducer  and  analog-to-digital  conversion  system 
required  to  accommodate  the  “upstream”  dynamic  power  range.  Consider  a  local  transient  such 
as  a  vacuum  pump  start,  which  has  a  5 A  current  amplitude.  Using  Equation  2-16  and  the  NILM 
configurations  in  TABLE  2-2,  the  effective  number  of  bits  Be,d  of  the  raw  current  for  the 
downstream  configuration  is 


Be,d  —  1°§2 


5  * 


1000 


*  155  *  212 


10 


2-19 


Be,d  ~  8.31  bits 


2-20 


Similarly,  the  effective  number  of  bits  Be,u  of  the  raw  current  in  the  upstream  is 

1 


Be,U  —  1°§2 


5  * 


2000 


*  5  *  212 


10 


2-21 


Be,u  ~  2-36  bits 


2-22 


This  reveals  approximately  6  bits  of  data  accuracy  is  lost,  which  results  in  significant 
degradation  of  the  raw  current  measured  signal.  Figure  2-8  shows  raw  data  measurements  for  a 
portion  of  a  vacuum  pump’s  transient  (shutoff)  taken  by  the  two  configurations  described  in 
TABLE  2-2.  Note  how  small  changes  of  current  in  the  aggregate  may  be  lost  in  the  raw  data 
measurements  due  to  the  quantization  that  occurs  due  to  reduced  accuracy. 
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Figure  2-8:  Two  NILM  configurations  recording  the  same  transient.  Top  plot  shows  raw  data  captured  on  both 
configurations  on  same  scale.  Bottom  plot  is  a  close-up  view  of  the  upstream  configuration  raw  data  capture  revealing 

quantization  error. 


2.2.2  Spectral  Envelope  Accuracy 

Current  NILM  analysis  methods  rely  solely  on  preprocessor  output,  not  raw  data  measurements. 
Proper  (5)  and  Branch  (7)  each  utilized  “prep”  output  to  perform  detection,  classification,  and 
diagnostic  analysis  for  both  the  CHT  system  and  Reverse  Osmosis  (RO).  More  specifically,  the 
first  harmonic  indices  of  real  and  reactive  power  have  proven  to  be  fruitful  signatures  for 
shipboard  NILM  applications. 
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Recall  from  Equations  2-17  and  2-18  that  Pi  and  Qi  are  calculated  using  the  time-average  over  a 
60  Hz  period.  Similar  to  a  low-pass  filter,  as  the  preprocessor  averages  the  raw  data  it  removes 
some  of  the  higher  frequency  content.  This  results  in  a  “smoother”  power  data  stream,  which 
consequently  reduces  the  raw  data  quantization  error.  Although  this  is  not  rigorously  proven, 
shipboard  NILM  aggregate  results  show  strong  evidence  of  this  “smoothing”  effect.  Based  on 
these  indications,  the  “prep  “output  is  less  susceptible  to  the  data  quantization  and  scaling  issues 
that  reduce  raw  data  accuracy. 

Figure  2-9  shows  the  preprocessor  output  comparison  between  the  downstream  and  upstream 
NILM  using  the  same  raw  data  compared  previously  in  Figure  2-8.  Although  the  raw  data  shows 
significant  quantization  in  the  upstream,  the  preprocessor  yields  very  similar  results  in  the 
upstream  and  downstream  output.  These  results  indicate  that  the  preprocessor  output  maintains 
the  power  information  necessary  to  perform  analyses  such  as  load  detection  and  diagnostics. 


Figure  2-9:  Preprocessor  spectral  envelope  comparison  for  raw  data  captured  in  Figure  2-8. 
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2.3  Measurement  Distortion 


Alteration  of  the  measured  voltage  and  current  waveforms  in  any  matter  may  adversely  impact 
NILM  performance.  Distortion  may  be  generated  internally  from  the  NILM  system  or  coupled 
from  external  sources;  this  falsification  of  measured  data  may  be  carried  throughout  the  data 
acquisition  and  processing  stages,  resulting  in  improper  classification  and  diagnostics. 
Unwanted  energy  such  as  thermally  generated  white  noise,  electromagnetic  coupling,  or 
magnetic  interference  must  be  mitigated  throughout  the  design  process. 

In  (1),  Paris  describes  two  general  categories  of  distortion  related  to  NILM  applications:  (1) 
uncorrelated  noise  and  (2)  correlated  noise  and  coupled  signals.  A  brief  description  of  these  two 
sources  is  described  for  the  reader;  however,  for  additional  information  on  these  topics  refer  to 
(1)  and  (16). 

2.3.1  Uncorrelated  Noise 

One  form  of  uncorrelated  noise  that  may  occur  within  the  system  and  which  serves  as  a  model 
for  additive  noise  in  general  is  additive  white  Gaussian  noise  (AWGN).  This  type  of  noise  is 
characterized  by  a  constant  spectral  density,  which  implies  the  energy  associated  with  the  noise 
is  evenly  distributed  across  all  frequencies.  Furthermore,  the  noise  is  normally  distributed 
around  a  mean  of  zero.  The  NILM  utilizes  a  Kalman  filter  to  reduce  error  associated  with  the 
zero-crossing  of  the  voltage  measurement.  The  robust  design  of  the  Kalman  filter  results  in  high 
noise  rejection  (15),  and  so  the  NILM  is  largely  unaffected  by  noise  in  the  voltage  input. 

Current  measured  by  the  preprocessor  is  also  susceptible  to  this  particular  noise.  However,  the 
NILM’s  spectral  envelope  calculations  remove  the  higher  frequency  content  of  the  incoming 
waveforms.  Since  the  noise  energy  is  equally  distributed  across  all  frequencies,  the  spectral 
output  only  contains  a  small  portion  of  the  energy  associated  with  the  noise.  Based  on  these 
informal  discussions  and  analysis  of  field  results,  it  is  highly  unlikely  that  this  type  of  noise  will 
impact  the  scalability  of  the  NILM. 
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2.3.2  Electromagnetic  Interference 


Noise  sources  external  to  the  NILM  may  also  impose  a  risk  to  its  overall  performance.  One 
potential  source  is  from  electromagnetically  coupled  sources.  Electromagnetic  energy  may  be 
generated  from  such  sources  as  current-carrying  wire  or  power  generators.  As  the  number  of 
loads  for  a  NILM  application  increases,  the  NILM’s  susceptibility  to  this  type  of  energy 
increases. 

One  method  for  mitigating  risk  associated  with  externally  coupled  energy  is  proper  shielding 
techniques.  The  NILM’s  enclosure  box  helps  isolate  the  data  acquisition  components  from  these 
externally-generated  emissions.  Additional  ways  to  reduce  susceptibility  include  smart  layout 
designs.  The  NILM’s  modular  design,  which  allows  physical  separation  between  data 
acquisition  components  and  processing  hardware,  is  also  an  important  aspect  in  electromagnetic 
compatibility.  Based  on  previous  field  studies,  which  includes  several  installations  on  Coast 
Guard  ships,  as  well  as  installations  at  the  LBES  facility,  the  NILM  does  not  appear  highly 
susceptible  to  electromagnetic  interference  in  practice. 

2.4  Disaggregation 

Perhaps  the  most  challenging  issue  facing  aggregate  power  monitoring  is  disaggregation  of  loads 
of  interest.  The  obvious  benefit  created  by  monitoring  more  loads  from  a  single  NILM  is  sensor 
count  reduction.  However,  this  imposes  new  challenges  when  it  comes  to  extracting  target  loads. 
The  complexity  associated  with  disaggregation  is  not  constant  from  one  site  installation  to  the 
next.  It  is  very  dependent  upon  the  type  of  information  that  must  be  obtained  in  order  to  fulfill 
the  requirements  for  a  particular  NILM  installation.  For  example,  one  NILM  application  may 
only  require  properly  identifying  the  status  of  multiple  pumps  within  a  fluid  system.  This 
scenario  could  most  likely  be  addressed  using  steady  state  power  levels  and  finite  state  modeling. 
A  different  application  may  require  loading  information  of  a  particular  system  such  as  the  low 
pressure  air  compressor  (LPAC)  system.  Since  LPAC  loading  results  in  variation  in  the 
compressor  cycling,  a  different  approach  than  that  used  for  pump  status  will  be  required. 
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2.4.1  Methods 


Disaggregation  is  the  act  of  separating  a  known  entity  into  component  parts.  In  order  for  the 
NILM  to  disaggregate  load  transients  within  the  collective  signal,  it  must  identify  a  unique 
feature  or  metric  that  distinguishes  the  known  entity  from  all  other  loads  and  background 
variations  within  the  signal.  Common  methods  used  to  detect  and  classify  individual  events  rely 
on  steady-state  power  changes,  finite  state  modeling,  and  transient  behavior.  These  approaches 
typically  utilize  first  order  power  and  reactive  power  to  achieve  the  desired  results.  More 
advanced  approaches  may  include  higher  order  harmonics,  as  well  as  frequency  analysis. 
Current  research  is  also  exploring  the  use  of  machine  learning  algorithms  that  are  capable  of 
analyzing  multiple  features  (i.e.  period,  amplitude,  frequency  content)  to  detect  specific  loads 
within  a  system. 

Recent  field  studies  have  relied  heavily  on  steady-state  power  levels  in  conjunction  with 
transient  behaviors  to  perform  detection  and  classification  of  loads.  Proper  (5)  used  these 
methods  to  perform  automated  classification  on  the  vacuum-assisted  sewage  system,  which 
includes  two  different  types  of  pumps  within  the  system.  Results  from  these  studies  showed  the 
validity  in  this  type  of  analysis.  As  the  NILM  scales  to  monitor  more  loads,  the  chance  for 
transient  overlap  increases,  which  could  adversely  affect  the  NILM’s  ability  to  extract  individual 
loads.  Furthermore,  higher  levels  of  system  monitoring  result  in  higher  overall  power  levels; 
these  levels  will  introduce  greater  levels  in  power  variation.  Ultimately,  the  increased  variations 
could  mask  the  transients  of  interest. 

Depending  upon  the  size  of  the  site  being  monitored,  the  overlap  issue  may  be  inconsequential. 
The  number  of  loads  may  be  small  enough  that  the  risk  associated  with  transient  overlap  may  be 
ignored.  For  larger  site  installations,  smarter  software  solutions  that  take  into  account  sequential 
steps  for  a  given  process  may  be  implemented.  For  instance,  Bennett  (4)  identified  the  NILM’s 
ability  to  identify  multiple  steps  within  a  gas  turbine  motor  start-up  procedure.  One  step  in 
particular  included  a  +28  VDC  signal  that  opened  a  vent  valve  downstream  of  the  exhaust  fan  for 
the  propulsion  turbine.  Although  the  actuator  signal  is  too  small  to  see  in  the  aggregate  level,  the 
associated  start-up  of  the  vent  fan,  which  occurs  nearly  immediately  following  the  actuator 
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signal,  could  be  used  to  classify  the  turbine  start-up.  In  general  cases  where  overlap  is  expected 
based  on  system  procedures,  the  multiple  events  may  be  treated  as  one  using  a  classifier  designed 
to  detect  the  multiple  events. 

Steady-state  power  variations  may  also  managed  using  techniques  similar  to  those  used  for 
overlap  issues.  Jones  (6)  identified  high  frequency  oscillations  associated  with  lube  oil  heaters 
during  his  research.  The  magnitude  of  these  oscillations  in  the  aggregate  is  capable  of  masking 
smaller  loads  of  interest.  One  solution  to  this  problem  is  to  use  a  filter,  such  as  a  notch  filter,  to 
remove  the  cycling  heater  masking  effect.  Another  useful  type  of  signal  processing  filter  is  the 
median  filter,  because  it  preserves  the  sharp  edges  of  the  power  signal,  which  is  often  used  for 
transient  identification. 

2.4.2  ESCANABA  Experimental  Data 

To  demonstrate  the  capability  of  the  NILM  to  disaggregate  loads,  two  monitoring  systems  were 
installed  on  the  USCGC  ESCANABA  in  September  2008.  An  “upstream”  NILM  (Figure  2-10) 
was  placed  at  P-2-282-3;  a  “downstream”  NILM  was  installed  at  the  CHT  control  panel.3  Data 
was  collected  over  a  two  month  period  while  the  ship  was  deployed,  then  retrieved  for  analysis 
in  November  2008. 


3  Refer  to  Section  2.2,  Measurement  Accuracy,  for  further  details  on  P-2-282-3. 
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Figure  2-10:  Upstream  NILM  System.  The  sensor  box  receives  the  raw  current  and  voltage  signals  from  the  panel, 
samples  the  data  at  8  kHz,  and  sends  the  resulting  data  to  the  computer  via  Ethernet. 


Based  on  the  data  acquisition  components  listed  in  TABLE  2-3  ,  the  dynamic  current  range  for 
the  upstream  NILM  is  larger  than  the  downstream  NILM  by  a  factor  of  2.  This  results  in  an 
effective  bit  increase  (Equation  2-16)  of  approximately  1  bit  from  upstream  to  downstream, 
which  does  not  introduce  significant  quantization  error. 


TABLE  2-3:  SENSOR  CONFIGURATION  FOR  DISAGGREGATION  EXPERIMENT 


Upstream 

Downstream 

Transducer 

LA-205 

LA- 150 

kn 

1/2000 

1/2000 

Rm 

50  Q 

100  Q 

The  objective  of  this  experiment  was  to  determine  whether  the  relatively  small  transients 
associated  with  the  CHT  vacuum  pump  could  be  reliably  identified  in  the  upstream  data.  Figures 
Figure  2-llthroughFigure  2-13  show  data  captured  from  October  22,  2008  from  0430-0500.  The 
downstream  plots  show  three  distinct  vacuum  pump  runs  with  no  additional  loads.  Upstream 
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plots  display  the  vacuum  pump  runs  along  with  additional  loads  on  a  steady  35  kW  power  level. 
The  closer  view  of  the  first  CHT  pump  run  (Figure  2-12)  displays  the  additional  transient  event 
in  the  upstream  during  the  CHT  pump  run.  Additionally,  it  illustrates  how  the  large  amplitude 
steady- state  variations  relative  to  downstream  data  may  prevent  proper  classification. 


0  5  10  15  20  25  .30 

Time  (min) 


Figure  2-11:  Full  window  time-aligned  comparison  of  downstream  data  (bottom)  and  upstream  (data).  Three  vacuum 
pump  transients  are  shown  in  the  downstream  (bottom).  The  three  vacuum  pump  transients  along  with  additional 
transients  and  noise  are  visible  upstream  (top)  on  top  of  a  steady  35  kW  load. 
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Figure  2-12:  Closer  view  of  the  first  CHT  pump  run,  with  a  larger  unrelated  transient  captured  in  the  upstream  data. 


Figure  2-13  illustrates  how  the  general  shape  of  the  downstream  transient  is  maintained  in  the 
aggregate.  Assuming  the  transient  shape  is  unique,  classification  may  still  be  possible  in  the 
upstream.  The  following  section  will  discuss  potential  solutions  to  load  classification  using 
transient  exemplars  or  “fingerprints”  to  match  with  similar  waveforms  in  the  aggregate. 
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Figure  2-13:  Detailed  view  of  the  start-up  transient  of  the  CHT  pump.  Note  the  general  shape  is  still  present  in  the 

upstream  data. 


2.4.3  Matching  Waveforms 


One  common  method  used  to  determine  the  similarity  in  two  sampled  waveforms,  /  [n]  and  g[n\, 
is  to  determine  the  Euclidean  distance,  D,  defined  as 

D2  =  ^  (f[n]  ~  g[n])2  2'23 

ueN 

When  expanded, 


D2 


^((/W)2  -  2 f[n]g[n]  +  ( g[n ])2) 
^/[n]2+  ^gin]2-  2^f[n]g[n] 


2-24 

2-25 


2-26 


As  an  inner  product,  this  becomes 

D2  =  (/•/)  +  (g-g)~  2(f-g) 
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2-27 


=  l/l2+  \g\2~  2(/  •  g) 


If  the  waveforms  match,  D  =  0,  and  the  equation  reduces  to 


0=  |/|2  +  \g\2~  2  {f.g) 


If  the  amplitudes  match,  |/|  =  \g\,  so 


f  •  g 


2  Iff  | 2 

2 


/  *  ff 

M2 


=  1 


2-28 


2-29 

2-30 


So,  given  an  input  signal,  /  [n],  and  an  exemplar,  g[n\,  if  the  two  waveforms  match  the  figure  of 
merit,  M,  approaches  1.0.  As  M  approaches  1,  it  can  be  reasonably  assumed  that  the  two 
waveforms  are  similar  in  amplitude  and  shape.  This  is  an  important  attribute  for  disaggregation 
methods  in  the  presence  of  many  differently-sized  loads  that  may  share  similar  shapes,  which  is 
often  the  case  with  shipboard  systems  that  contain  multiple  induction  motors. 


When  an  exemplar  is  chosen  for  a  given  load,  it  will  contain  a  minimal  number  of  indices  to 
capture  the  load’s  characteristic  of  interest.  This  exemplar  can  “slide”  across  the  data  stream  and 
a  value  for  M  computed.  As  the  exemplar  approaches  a  feature  in  the  data  stream  that  matches 
in  shape  and  amplitude,  the  metric  M  will  approach  1.0. 


A  faster  method  of  computing  the  same  metric  is  based  on  cross-correlation  (16).  Given  two 
waveforms  x  and  y,  cross-correlation  is  defined  as: 


(x  *  y)  [n] 


OO 

^  x*  [m]y[n  +  m] 

m=  -oo 


2-31 


where  x*  denotes  the  complex  conjugate  of  x.  For  real  numbers,  this  represents  the  same  basic 
operation  as  the  dot  product  (x-y).  Cross-correlation  can  be  computed  using  the  discrete 
Fourier  transform  F  using  the  relation 

F{x*y}=rrx}F{y}  2-32 
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Application  and  Results 


The  cross-correlation  method  was  attempted  using  data  from  the  patrol  period.  The  candidate 
transient  was  a  vacuum  pump  startup.  First,  an  exemplar  was  extracted  from  the  downstream 
data  sample  and  the  DC  offset  was  removed  (Figure  2-14). 


0  0.05  0.1  0.15  0.2  0.25  0.3 


Time  (s) 

Figure  2-14:  Exemplar  used  to  match  CHT  vacuum  pump  startup  transient. 

Then,  the  correlation  was  computed  using  the  exemplar  and  the  full  window  of  data  from  Figure 
2-11.  First,  the  downstream  data  was  analyzed.  Figure  2-15  reveals  a  positive  match  between 
the  exemplar  and  downstream  data.  As  expected,  the  correlation  value  approaches  1.0  as  the 
exemplar  “slides”  over  the  transient  in  the  downstream  data  stream.  A  peak-detection  algorithm 
with  a  threshold  of  ±  0. 1  was  used  in  this  application. 
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Figure  2-15:  Positive  exemplar  match  in  the  downstream  data. 
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Next,  the  exemplar  extracted  from  the  downstream  data  was  used  to  perform  cross-correlation  in 
the  upstream  data.  Figure  2-16  illustrates  the  results.  Again,  the  method  correctly  identified  the 
transient  of  interest  even  with  the  steady-state  35  kW  power  level  and  other  transient  events. 


<D 

D 

? 

c 

o 


<L> 

g 

o 

U 


2 

1.5 

1 

0.5 

0 

-0.5 

-1 


2.78 


2.79 


2.82 


Upstream  correlation  - 

Detection  window  -  - 

- 

- 

_  i 

1 

2.83 


2.80  2.81 
Time  (min) 

Figure  2-16:  Positive  exemplar  match  against  a  CHT  pump  startup  in  the  upstream  data. 


Finally,  Figure  2-17  shows  an  example  of  the  correlation  method  correctly  distinguishing  a  non- 
CHT  pump  startup  transient  in  the  upstream.  Although  the  correlation  value  approaches  and 
passes  through  the  peak  detection  window,  the  peak  occurs  at  a  higher  correlation  value,  which  is 
easily  distinguishable  using  the  peak  detection  algorithm. 


Figure  2-17:  Positive  exemplar  mismatch  against  a  non-CHT  event. 
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For  this  particular  case,  the  cross-correlation  technique  correctly  identified  all  pump  matches. 
Similar  analysis  was  conducted  on  21  disjoint  hours  over  a  three-week  period  from  October  10  to 
October  20.  For  every  CHT  pump  transient  found  in  the  clean  downstream  data,  the  same  event 
was  correctly  identified  in  the  upstream  data.  Furthermore,  there  were  no  false  positives  or 
negatives,  indicating  that  this  method  can  quickly  and  reliably  disaggregate  target  loads  even  in 
the  presence  of  large  variations  and  unrelated  transients. 
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Chapter  3 


LBES  Aggregate  Power  Analysis 


Methods  of  disaggregation  are  examined  using  field  data  collected  from  the  Navy’s  Land  Based 
Engineering  Site  (LBES)  located  in  Philadelphia,  Pennsylvania.  The  facility  contains  a  fully 
functional  ARLEIGH  BURKE  (DDG-51)  main  machinery  space,  which  is  used  by  the  navy  for 
crew  certifications  and  also  serves  as  a  test  site  for  various  mechanical,  electrical,  and  software 
upgrades.  NILM  research  has  been  conducted  at  the  LBES  site  over  the  past  few  years  with 
emphasis  on  local  diagnostic  indicators  and  supervisory  control  ( (17)  (6)  (4)  (18)). 

LBES  studies  began  with  NILM  installations  on  the  fuel  oil  (FO)  and  low  pressure  air 
compression  (LPAC)  systems  in  2005  (18).  This  was  followed  with  additional  NILM 
installations  to  monitor  the  General  Electric  LM  2500  gas  turbines  used  for  propulsion,  along 
with  various  supervisory  control  systems  associated  with  the  turbines  (17).  Additional  loads 
monitored  through  NILM  at  LBES  include  fuel  oil  service  pumps  (FOSPs),  lube  oil  service 
pumps  (LOSPs),  and  gas  turbine  generators  (GTGs).  In  2007,  Jones  (6)  installed  two  NILM 
systems  on  the  main  switchboards  supplying  power  to  all  of  the  auxiliary  system  loads 
previously  mentioned.  In  June  2008,  data  was  collected  from  the  switchboard  NILMs  and 
various  service  loads.  The  focus  of  this  section  is  to  examine  how  well  cross  correlation  and 
peak  detection  methods  perform  disaggregation  of  the  various  downstream  loads  in  the  upstream 
data. 

3.1  Aggregate  (3SA)  Description 

Although  two  NILMs  were  installed  upstream  at  the  LBES  site,  only  one  (3SA)  will  be  used  for 
the  purpose  of  this  study,  since  all  target  loads  are  on  3SA.  The  CT  chosen  for  this  installation 
differs  from  the  previously  described  CTs  (i.e.  LA  55-P,  205-S).  This  particular  CT,  the  Fluke 
i3000,  is  a  split-core  CT  that  clamps  around  the  current-carrying  wire.  The  decision  to  use  the 
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clamp  style  was  made  based  on  the  site  engineer’s  request  to  avoid  breaking  electrical 
connections  in  the  switchboard  (6) . 

Selection  of  the  Fluke  i3000  resulted  in  the  necessity  to  use  an  amplifier  circuit  to  increase  the 
signal  level  into  the  PCI-1710  ADC.  The  CT  offers  three  current  ranges:  30A,  300A,  and 
3000A.  Based  on  potential  currents  at  the  switchboard,  the  3000A  range  was  chosen.  Unlike 
typical  NILM  CTs,  this  particular  device  converts  the  primary  current  IP  to  voltage  using  KN  = 
0.1  mY/A;  there  is  no  burdened  resistor  required  upstream  of  the  ADC.  This  results  in  extremely 
small  output  voltages  that  require  amplification.  Amplification  was  achieved  using  an  amplifier 
card  designed  by  Professor  Robert  Cox  (6).  The  amplifier  card  receives  Is  and  splits  the  signal 
into  two  levels,  Channels  1  and  2.  The  output  of  channel  1  is  amplified  by  a  factor  of  10,  while 
the  output  of  channel  2  is  amplified  by  a  factor  of  100.  For  the  purposes  of  upstream  NILM 
discussion  at  LBES,  the  two  outputs  of  the  amplifier  cards  will  be  called  measured  voltage 
ranges  Vm .  This  follows  data  acquisition  nomenclature  defined  in  previous  sections  for  NILMs 
equipped  with  burdened  resistors. 

A  sample  power  plot  is  provided  in  Figure  3-1  to  illustrate  the  relative  difference  in  scaling 
between  channels  1  and  channel  2  data.  There  are  two  factors  that  influence  the  scale  difference. 
First,  the  obvious,  is  the  difference  in  amplification  of  each  signal  prior  to  reaching  the  ADC. 
Second,  the  designers  chose  to  select  different  gains  for  each  channel.  Recall  the  PCI- 17 10  card 
was  used  for  this  installation.  This  card  offers  selectable  voltage  ranges,  Vr.  Channel  1  was  set 
to  +/-  2.5  V;  channel  2  was  set  to  +/-  10V.  Based  on  these  scale  factors,  raw  count  output  from 
the  ADC  of  channel  2  is  approximately  2.5  times  that  from  channel  lfor  a  given  IP.  A  summary 
table  of  the  two  NILM  arrangements  is  provided  in  TABLE  3-lbelow.  This  chapter  will  use  the 
channel  1  data  only,  since  it  has  a  higher  degree  of  scaling  involved,  which  yields  the  most 
conservative  analyses. 
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Time  (  Seconds  ) 

Figure  3-1:  Power  plot  of  both  channels  available  for  3SA.  The  direct  current  (d.c.)  bias  was  removed  from  each  signal 
to  provide  a  comparative  illustration  of  the  scale  difference  between  the  two  channels. 


TABLE  3-1:  3SA  NILM  CONFIGURATION 


Upstream  CH.l 

Upstream  CH.2 

Current  Transducer 

Fluke  i3000 

Fluke  i3000 

kn 

.1  mV/A 

.1  mV/A 

Rm 

N/A 

N/A 

G 

10 

100 

VR 

+/-  2.5  V 

+/-  10  V 

Bits 

12 

12 

3.2  L0SP2A 

Two  LOSPs,  LOSP  2A  and  2B,  are  piped  to  a  common  header  supply  lubricating  oil  the  main 
reduction  gears  and  the  shaft  thrust  bearing.  Each  pump  is  a  two-speed  variety  with  motor 
ratings  of  30/60  HP  (slow/fast).  LOSP  2A  is  monitored  using  a  200  burdened  resistor  and  a 
LEM  LA-305S  current  transducer.  The  LA-305S  has  a  KN  value  of  1/2500  with  a  maximum 
primary  current  range  of  +/-  300  amps.  The  NILM  uses  a  PCI- 17 10  ADC  with  a  VR  of  +/-  5V. 
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TABLE  3-2:  LOSP  2A  NILM  CONFIGURATION 


LOSP  2A 

Current  Transducer 

LEM  LA-305 S 

KN 

1/2500 

Rm 

20  Q 

Vr 

+/-  5V 

Bits 

12 

Recall  from  Chapter  2  that  the  effective  bits  calculation  provides  a  good 
quantization  from  the  local  NILM  to  further  upstream.  The  effective  bits, 
3SA,  BE,um&y  be  written  as 


Be,u  —  1°§2 


I  *  KN  *  G  *  2 
Wr\ 


BN 


measure  of  the 
Be,  for  channel 


data 
1  of 


3-1 


where  G  is  the  gain  from  the  amplifier  card.  Using  an  IP  of  +/-  72  amps  for  the  LOSP4,  BF  u  for 
3SA  channel  1  is 


Be,u  —  1°§2  I  " 


(12  *  le  —  04  *  10  *  212N 


\ 


3-2 


»  5.882  bits 


3-3 


Equation  2-16  yields  the  following  downstream  Be,d- 

(72  *  1/2500  *  20  *  212N 


Be,d  —  1°§2  I  " 


\ 


10 


3-4 


«  7.965  bits 


3-5 


The  difference  in  effective  bits  between  the  channel  1  data  Be,u  and  the  downstream  LOSP  2 A 
data  Be, pis  only  2.2  bits. 


4  Full  load  amperage  rating  for  high  speed  operation  (4). 
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3.2.1  LOSP  2A  OFF  to  SLOW  Cross  Correlation  Results 


Data  obtained  on  June  19,  2008,  between  1030  and  1200  was  analyzed.  During  this  time  period, 
four  locally-monitored  loads  were  cycled  along  with  additional  transients  in  the  aggregate. 
Figures  3-2  through  3-4  provide  comparisons  of  the  time-aligned  spectral  power  envelopes 
captured  both  upstream  and  downstream.  Cross  correlation  and  peak  detection  methods 
described  in  Chapter  2  were  used  to  extract  the  LOSP  2A  OFF/SLOW  transient  from  the 
aggregate. 


Figure  3-2:  Time-aligned  power  spectral  plots  of  3SA  (upstream)  and  LOSP  2A  (downstream).  Note  the  aggregate 
contains  additional  load  starts  including  LOSP  2B  transients,  which  are  highlighted  in  the  follow-on  figure. 
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Figure  3-3:  Closer  view  of  aggregate  and  downstream  power  plots.  The  aggregate  contains  LOSP  2A  and  LOSP  2B 
transients.  The  lower  plot  shows  LOSP  2A  cycle  through  its  various  states  of  operation  (OFF/SLOW/FAST). 


Figure  3-4:  Detailed  view  of  the  startup  transient  from  OFF  to  SLOW  speed.  Aggregate  contains  a  higher  frequency 
cycling  load  that  was  identified  to  be  Lube  Oil  heaters  by  Jones  (6). 

In  order  to  conduct  the  cross  correlation  technique,  a  detailed  exemplar  of  the  transient  of  interest 
is  required.  Figure  3-4  shows  a  detailed  view  of  the  downstream  transient  of  interest  along  with 
the  time-aligned  upstream  data.  The  sample  exemplar  for  OFF  to  SLOW  Speed  was  extracted 
from  the  downstream  data,  and  then  the  DC  offset  was  removed.  The  resulting  template  is 
shown  in  Figure  3-5. 
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Figure  3-5:  Exemplar  used  for  LOSP  2 A  OFF  to  SLOW  transient. 


Next,  cross-correlation  was  computed  using  this  exemplar  against  full  hour  windows  of 
downstream  and  upstream  data.  Recall  from  Chapter  Two,  that  a  strong  correlation  is  indicated 
by  a  local  maxima  near  y=  1.0  for  a  given  time  t.  A  peak  detection  window  between  0.9  and  1.1 
was  used  as  the  threshold  for  a  positive  match.  Figures  3-6  and  3-7  show  the  results  for  the 
downstream  correlation.  Each  OFF  to  SLOW  event  was  correctly  identified  and  all  other  events 
fell  outside  of  the  detection  range. 
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Figure  3-6:  Detection  window  of  downstream  correlation  with  multiple  transient.  The  two  spikes  within  the  threshold 
values  were  both  correctly  identified  LOSP  starts  (OFF  to  SLOW).  The  spikes  below  the  threshold  window  correspond  to 

SLOW  to  FAST  and  FAST  to  OFF  events. 


Figure  3-7:  Exemplar  match  against  a  LOSP  start  transient  recorded  in  the  downstream. 


After  downstream  analysis,  the  exemplar  was  cross  correlated  in  the  upstream  data.  Figure  3-8 
shows  a  positive  match  for  the  LOSP  2 A  OFF/SLOW  transient  in  the  aggregate.  Figure  3-9 
shows  a  positive  mismatch  that  corresponds  to  a  LOSP  2A  SLOW/FAST  transient.  All  LOSP 
2A  OFF/SLOW  transients  were  correctly  identified  in  the  aggregate  with  no  false  positives 
or  negatives. 
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Figure  3-8:  Upstream  correlation  window  for  positive  LOSP  2 A  OFF/SLOW  match. 


Figure  3-9:  Positive  mismatch  in  upstream  data  for  LOSP  2A  OFF/SLOW. 


3.2.2  LOSP  2A  SLOW  to  FAST  Cross  Correlation  Results 


The  transient  due  to  a  shift  in  speed  from  SLOW  to  FAST  of  the  LOSP  was  scrutinized  next.  An 
exemplar  was  extracted  from  the  downstream  data.  The  same  data  samples  used  for  the  previous 
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discussion  on  OFF  to  SLOW  were  used  for  this  study.  Refer  to  Figure  3-2  for  the  power  plot 
results.  The  exemplar  chosen  for  this  transient  is  provided  as  Figure  3-10. 


Figure  3-10:  Exemplar  used  to  match  SLOW  to  FAST  transient  of  LOSP  2A. 


Figures  3-11  and  3-12  are  results  from  the  cross-correlation  method  using  downstream  data  and 
upstream  data,  respectively.  Again,  the  detection  method  correctly  identified  the  transient  of 
interest  both  downstream  and  upstream.  There  were  no  false  positives  or  negatives  revealed. 


Figure  3-11:  Positive  match  between  exemplar  and  downstream  LOSP  2A  SLOW/FAST  transient.  Each  candidate 
transient  (SLOW/FAST)  was  correctly  identified;  there  were  no  false  positives. 
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Figure  3-12:  Positive  match  between  exemplar  and  upstream  LOSP  2A  SLOW/FAST  transient. 


3.3  Low  Pressure  Air  Compressor  (LPAC)  Number  2 

The  LPAC  provides  compressed  air  at  approximately  125  psig  to  the  LBES  low  pressure  air 
system.  The  low  pressure  air  serves  to  operate  various  valves  and  provide  pneumatic  power  for 
certain  plant  equipment.  The  compressor  is  an  Ingersoll-Rand  NAXI  100-4A  driven  by  a  35  HP 
squirrel  cage  motor  rated  for  35  full  load  amps  (4).  NILM  summary  is  listed  in  TABLE  3-3. 
Using  the  35  amp  full  load  rating,  BE,u  =  4.841  bits  and  Be,d  =  7.689  bits. 
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TABLE  3-3:  NR  2  LPAC  NILM  CONFIGURATION5 


Nr  2  LPAC 

Current  Transducer 

LEM  LA-305 S 

KN 

1/2500 

Rm 

36  Q 

Vr 

+/-  5V 

Bits 

12 

3.3.1  LPAC  Nr  2  Correlation  Results 

Figure  3-13  is  a  time-aligned  spectral  power  plot  of  downstream  and  upstream  data  recorded  on 
June  19,  2008  between  1100  and  1200.  Two  LPAC  Nr  2  compressor  start-ups  and  shut-downs 
were  recorded  in  the  upstream  and  downstream  data  during  this  window.  This  performance  does 
not  represent  normal  operation  of  the  LPAC  unit,  but  was  conducted  to  provide  known  transients 
for  research.  During  normal  operation,  the  compressor  remains  energized  and  the  air  system  is 
under  load.  (6)  Figure  3-14  is  a  close-up  view  of  a  complete  compressor  cycle  captured  locally 
and  upstream. 


5  A  pressure  transducer  was  also  installed  to  provide  compressor  discharge  pressure  measurements.  However,  this  is  not  relevant 

to  this  discussion. 
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Figure  3-13:  Time-aligned  full  power  plot  of  3SA  and  local  LPAC  Nr  2  NILM.  Two  transients  on  bottom  plot  are 

compressor  cycles  from  OFF  to  ON,  then  back  OFF. 


Figure  3-14:  Closer  view  of  LPAC  Nr  2  compressor  cycle. 
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The  cross-correlation  method  was  used  on  the  downstream  data,  followed  by  the  upstream  data. 
The  exemplar  in  Figure  3-15  was  extracted  from  the  downstream  data.  The  downstream  data 
yielded  positive  exemplar  matches  for  all  target  transients  (Figure  3-16);  however,  there  were 
many  false  positives  identified  in  the  aggregate  (Figure  3-17  and  Figure  3-18).  Several  different 
exemplars  were  extracted  from  the  downstream  data  in  an  attempt  to  find  a  better  “fingerprint”. 
However,  the  “cleanest”  results  for  the  upstream  aggregate  were  derived  using  the  exemplar  in 
Figure  3-15. 


Figure  3-15:  LPAC  Nr2  exemplar  extracted  from  downstream  data. 
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Figure  3-16:  Positive  exemplar  matches  with  two  LPAC  compressor  start-ups  in  the  downstream  data. 
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Figure  3-17:  Close-up  view  of  the  upstream  correlation  results.  Data  captured  around  the  actual  target  transients. 


58 


1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 

- Upstream  Correlation 

- Detection  Window 


Time  (min) 

Figure  3-18:  Detailed  view  of  upstream  correlation  results.  All  correlation  values  within  detection  window  were  falsely 

identified  as  positives. 
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3.4  Fuel  Oil  Service  Pump  2A 


The  fuel  oil  service  pumps  (FOSPs)  supply  fuel  to  the  gas  turbine  motors  and  to  the  gravity  drain 
tank  used  by  the  gas  turbine  generators.  Each  FOSP  is  coupled  to  a  two  speed  motor.  During 
normal  operations,  only  one  pump  is  required  to  meet  system  demand.  FOSP  2A  is  monitored 
locally  by  a  NILM  with  configurations  listed  in  TABLE  3-4.  Based  on  a  10  A  full  load  rating  for 
the  FOSP  motor,  the  Be,d=  7.965  bits  and  BE,u=  3.034  bits,  which  is  nearly  5  bits  difference. 
This  would  result  in  significant  quantization  in  the  raw  data;  however,  preliminary  studies  in 
Chapter  2  indicate  the  preprocessor  output  is  largely  unaffected  by  the  quantization. 


TABLE  3-4:  FOSP  2A  NILM  CONFIGURATION 


FOSP  2A 

Current  Transducer 

LEM  LA-55S 

KN 

1/1000 

Rm 

61  Q 

Vr 

+/-  5V 

Bits 

12 

Similar  to  the  results  for  the  LPAC  Nr  2  correlation,  the  cross  correlation  method  correctly 
labeled  transients  in  the  downstream;  however,  it  failed  to  correctly  label  FOSP  2A  transients  in 
the  upstream.  Figures  3-19  and  3-20  provide  time-aligned  spectral  plots  of  both  downstream  and 
upstream  data.  It  seems  fairly  evident  from  the  close-up  view  that  the  FOSP  transient  cannot  be 
“seen”  by  the  naked  eye  due  to  the  large  cycling  loads  present  in  the  aggregate. 
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Figure  3-19:  Time  aligned  plots  of  upstream(top)  and  downstream(bottom)  data.  Transients  in  bottom  plot  are  FOSP  2 A 
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Figure  3-20:  Close-up  view  of  the  power  spectral  plots  for  downstream  and  upstream  data.  The  cycling  load  upstream 

masks  the  FOSP  transient. 


The  exemplar  chosen  for  the  FOSP  2A  OFF/SLOW  transient  is  provided  in  Figure  3-21. 
Various  exemplars  were  examined  to  include  leading  edge  only,  trailing  edge  only,  and  full 
exemplars.  Each  yielded  similar  poor  results  when  cross-correlated  with  the  upstream  data.  The 
downstream  cross-correlation  yielded  positive  matches  and  mismatches,  with  no  false 
identifications.  A  positive  match  is  shown  in  Figure  3-22. 
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Figure  3-21:  FOSP  2A  OFF/SLOW  exemplar  extracted  from  downstream  data. 


Figure  3-22:  Positive  cross-correlation  match  for  a  FOSP  2A  OFF/SLOW  transient. 


The  cross-correlation  results  for  the  upstream  are  illustrated  in  Figures  3-23  and  3-24.  Notice 
that  the  cross-correlation  value  does  approach  1.0  near  the  target  transient;  however,  several  false 
positives  are  also  generated. 
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Figure  3-23:  Close-up  view  of  cross-correlation  results  in  the  aggregate  data  stream.  The  FOSP  2A  OFF/SLOW  transient 

is  indicated. 
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Figure  3-24:  Cross-correlation  results  in  upstream  data  showing  false  positive  identifications  (highlighted  in  red). 


3.5  Gas  Turbine  Propulsion  Main  Engine  Controller 

The  gas  turbine  propulsion  main  engine  (GTM)  controller  NILM  monitors  a  total  of  four  loads. 
The  first  load,  power  panel  P  1-282-1,  receives  440  VAC  and  steps  it  down  to  120  VAC.  This 
panel  supplies  the  remaining  three  loads,  which  include  two  Universal  Engine  Controllers 

(UECs)  and  the  Shaft  Control  Unit.  The  UECs  primary  function  is  to  provide  control  and 
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monitoring  of  the  GTM  to  include  such  functions  as  power  level  control  and  overspeed 
protection.  The  SCU  primarily  provides  local  control  to  shaft  components,  the  controllable 
reversible  pitch  propeller,  and  other  auxiliary  equipment.  For  a  detailed  explanation  of  these 
loads  refer  to  (4).  The  NILM  is  arranged  with  one  voltage  tap  on  the  power  panel  and  four 
current  channels  for  each  load  (TABLE  3-5). 


TABLE  3-5:  GTM  NILM  CONFIGURATION 


Panel 

UEC  2A  &  2B 

SCU 

Current  Transducer 

LEM  LA-205 S 

LEM  LA-55S 

LEM  LA-55S 

kn 

1/2000 

1/1000 

1/1000 

Rm 

67  Q 

130  Q 

130  Q 

Vr 

+/-  2.5  V 

+/-  .625  V 

+/- 1.25  V 

Bits 

12 

12 

12 

This  discussion  will  focus  on  a  normal  GTM  start  sequence.  The  “upstream”  panel  will  be  the 
120  VAC  panel  P-1-282-2,  which  supplies  power  to  both  UECs  and  the  SCU.  Current 
amplitudes  associated  with  transients  generated  from  the  SCU  and  UECs  are  approximately  1  to 
3  amps.  Using  2A,  the  BE,u~  0.712  bits.  In  other  words,  the  signals  are  too  small  for  aggregate 
monitoring.  A  logical  alternative  from  a  design  perspective  is  to  use  the  power  panel  that 
directly  feeds  the  load.  Using  2A  for  /,  the  Be, u, panel  -  5.778  bits.  Be,d,uec  ~  9.734  bits  and 
BE, d, scu  ~  8.734  bits. 

A  normal  GTM  start  sequence  is  illustrated  in  Figure  3-25.  The  process  begins  with  the  opening 
of  the  starter  air  regulating  valve  (SARV),  which  is  accomplished  via  a  +  28  VDC  signal 
supplied  from  the  UEC.  Next,  the  UEC  opens  the  fuel  oil  service  valves  (FOSVs)  to  allow  fuel 
supply  to  the  engine.  At  the  same  time,  the  GTM  igniters  are  energized  (shown  in  Panel).  The 
SARV  is  closed  next  by  the  UEC.  The  SCU  supplies  a  +28  VDC  signal  to  a  relay  to  open  the 
GTM  vent  damper,  while  the  GTM  cooling  fan  is  energized.  The  cooling  fan  has  a  130  HP 
motor  and  produces  a  significant  voltage-sag,  which  is  seen  in  all  three  power  plots. 
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Figure  3-25:  Time-aligned  spectral  plots  of  a  normal  GTM  start  sequence.  The  large  sub- transient  voltage  sag  is  a  result 

of  the  GTM  cooling  fan  start.  (4) 


3.5.1  GTM  Cross-Correlation  Results 

SCU  Transient 

First,  the  SCU  was  chosen  as  the  downstream  candidate.  The  transient  generated  by  the  SCU  is 
due  to  the  +28  VDC  signal  sent  to  a  relay  to  open  the  GTM  vent  damper.  Recall  that  this  is 
immediately  followed  by  a  large  voltage  droop  due  to  the  start-up  of  the  GTM  cooling  fan.  The 
spectral  plots  for  the  upstream  and  downstream  data  are  provided  in  Figure  3-26.  Each  plot 
contains  significant  higher  frequency  content  whose  origin  is  unknown. 
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Figure  3-26:  Time-aligned  plots  of  Power  Panel  and  SCU.  Transient  of  interest  is  the  relay  signal  provided  by  the  SCU  to 

open  the  GTM  vent  damper. 


Several  iterations  of  exemplars  were  examined  using  methods  previously  described.  Ultimately, 
the  exemplar  shown  in  Figure  3-27  was  chosen.  Notice  the  footprint  contains  many  peaks  and 
valleys  and  does  not  provide  a  smooth  curve. 


Figure  3-27:  Exemplar  for  GTM  damper  valve  signal. 
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Downstream  and  upstream  cross-correlation  was  conducted  using  the  exemplar.  Results  are 
provided  in  Figures  3-28  and  3-29  below.  Downstream  yielded  positive  results;  however,  the 
upstream  failed  to  correctly  identify  the  event. 
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Figure  3-28:  Downstream  cross-correlation  results.  Even  with  significant  "chatter",  the  cross-correlation  method 

correctly  identified  the  transient. 


2 


Time  (min) 


Figure  3-29:  Upstream  cross-correlation  results.  Method  failed  to  correctly  identify  the  transient  and  generated  false 

positives. 
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Results  indicate  that  the  higher  frequency  “chatter”  severely  limits  the  cross-correlation  detection 
method.  A  simpler  approach  based  on  finite  state  modeling  and  steady  state  power  levels  would 
be  a  viable  alternative. 

UEC  SARV  Transient 


Transients  associated  with  the  UEC  are  captured  in  the  upstream  and  downstream  NILMs  in 
Figure  3-30.  The  first  transient  tested  is  the  relay  signal  provided  to  open  the  SARV.  This  signal 
only  produces  a  5  W  change  in  power,  which  is  the  smallest  transient  examined  thus  far. 
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Figure  3-30:  Time-aligned  plots  of  Panel  (aggregate)  and  UEC  (downstream). 


The  exemplar  for  the  SARV  signal  was  extracted  from  the  downstream  data  (Figure  3-31).  The 
downstream  data  yielded  positive  cross-correlation  results  (Figure  3-32).  The  disaggregation 
from  the  upstream  data  failed  again.  Figure  3-33  shows  how  the  correlation  between  the 
waveforms  produced  multiple  false  positives.  The  SARV  opening  transient  occurs  around  39.08 
minutes. 
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Figure  3-31:  Exemplar  of  UEC  SARV  transient  extracted  from  downstream  data. 
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Figure  3-32:  Downstream  cross-correlation  values.  The  transient  was  correctly  identified. 
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Figure  3-33:  Upstream  Correlation  results  for  UEC  SARV  exemplar.  Transient  of  interest  is  embedded  in  spikes  around 

39.08  minutes. 


UEC  FOSV  Transient 


Next,  the  FOSV  transient  (see  Figure  3-30)  was  extracted  from  the  downstream  data  and 
produced  the  exemplar  shown  in  Figure  3-34. 


Figure  3-34:  UEC  FOSV  exemplar  extracted  from  downstream  data. 
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Cross  correlation  results  are  shown  in  Figures  3-35  and  3-36.  Downstream  results  were 
satisfactory.  Upstream  results  revealed  strong  correlation  between  the  exemplar  and  the 
upstream  FOSV  transient.  However,  multiple  false  positives  were  captured. 


Figure  3-35:  Downstream  correlation  results  for  FOSV  transient.  Unique  waveform  and  amplitude  produced  strong 

correlation  relative  to  background  “chatter”. 
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Figure  3-36:  Upstream 


correlation  results  for  FOSV  transient.  Positive  match  for  event;  however,  there  were  multiple 


false  positives  as  well. 
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Chapter  4 


CHT  Classification  and  Diagnostics  Using  Ginzu  NILM 


The  first  operational  shipboard  NILM  application  was  achieved  by  Ramsey  in  2004.  (2)  This 
installation  was  followed  by  various  trend  analyses  of  cycling  fluid  systems  such  as  CHT  and 
ASW.  (19)  Further  NILM  development  focused  on  creating  robust  models  that  enhanced 
detection  of  improper  pump  operations  and  leaks  within  fluid  systems.  (3)  In  recent  years, 
research  teams  have  continued  to  gather  shipboard  system  data  in  efforts  to  improve  detection 
and  diagnostic  capabilities.  Initially,  this  data  was  stored  on  the  NILM  computer,  and  then 
processed  in  the  laboratory  to  determine  system  characteristics.  In  2008,  the  first  Ginzu  NILM6 
configuration  capable  of  real-time  detection,  classification,  and  system  diagnosis  was  installed. 

(5) 

In  September  2008,  the  USCGC  ESCANABA  deployed  with  the  Ginzu  NILM  installed  on  the 
CHT  system.  This  marks  the  second  patrol  period  where  the  Ginzu  NILM  was  installed.  This 
chapter  will  describe  the  key  characteristics  of  the  Ginzu  NILM  along  with  some  of  its 
capabilities.  The  methods  used  for  event  detection,  transient  classification,  and  transient 
diagnostics  are  also  discussed.  Data  captured  from  this  period  is  scrutinized  to  determine  the 
overall  effectiveness  of  the  NILM  with  regards  to  trend  analyses  and  diagnostics. 

4.1  CHT  System  Description 

The  CHT  system’s  function  is  to  collect  store,  and  transfer  waste  from  twenty  one  locations 
throughout  the  ship.  Central  to  the  system  is  the  360  gallon  vacuum  collection  tank  (VCT). 
Drains  from  the  ship’s  urinals,  toilets,  and  garbage  grinder  are  fed  to  the  VCT.  Two  alternative 
vacuum  pumps  energize  to  maintain  the  appropriate  vacuum  within  the  tank.  Pressure  switches 


6  The  Ginzu  NILM  was  the  NILM  configuration  used  by  Proper  (5)  in  January  2008  on  board  the  USCGC 
ESCANABA  to  monitor  the  RO  and  CHT  system. 
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located  atop  the  tank  are  used  to  provide  low  vacuum  signals  to  the  control  panel.  At  14  in-Hg,  a 
single  vacuum  pump  is  operated  to  evacuate  the  air  from  the  VCT.  If  vacuum  reaches  12  in-Hg, 
both  pumps  are  energized  to  increase  vacuum  level  within  the  tank.  When  the  pressure  switches 
sense  18  in-Hg,  all  operating  vacuum  pumps  are  secured.  Two  alternating  discharge  pumps  are 
used  to  remove  waste  collection  from  the  collection  tank  to  a  holding  tank  when  a  high  level 
setpoint  probe  is  energized;  the  pumps  are  secured  at  a  low  level  limit.  The  local  control  panel 
receives  its  inputs  from  the  various  sensors  and  actuates  the  appropriate  pumps  to  restore  normal 
operating  parameters.  The  NILM  monitors  all  four  pumps  via  the  local  control  panel  power 
signals.  A  basic  system  diagram  is  illustrated  in  Figure  4-lbelow. 


Figure  4-1:  NILM  Block  Diagram.  Dashed  lines  indicate  electrical  connections  within  the  system.  Solid  lines  represent 

fluids  (air  and  sewage).  (20) 


4.2  Ginzu  NILM  Description 

Figure  4-2  illustrates  the  core  functions  of  the  hardware  and  software  associated  with  the  Ginzu 
NILM.  Raw  current  and  voltage  measurements  are  collected  at  a  sampling  rate  of  8  kHz  via 
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voltage  and  current  sensing  lines.  This  raw  continuous  data  is  collected  by  the  NILM  sensor  box, 
where  it  is  converted  into  a  quantized  digital  signal  as  it  passes  through  an  ADC  card. 


The  quantized  digital  signal  from  the  ADC  is  passed  via  Ethernet  to  the  NILM  computer  for 
follow-on  processing.  This  processing  includes  conversion  of  the  high  frequency  signal  into 
power  envelopes  (discussed  in  Chapter  2),  as  well  as  detection  and  classification  techniques  used 
to  indicate  system  status  and  health  information.  At  the  heart  of  the  Ginzu  NILM  is  the  software 
package,  Ginzu.  The  basic  program  flow  consists  of  detection,  classification,  and  verification. 
In  its  normal  mode  of  operation,  Ginzu  receives  preprocessor  data  at  a  frequency  of  120  Hz.  Ten 
second  data  windows  are  loaded  and  passed  to  an  event  detector,  which  identifies  rapid  power 
changes.  If  an  event  has  occurred,  the  classifier  algorithm  is  used  to  identify  the  transient  based 
on  power  levels  and  system  state  information.  The  event  files  created  by  Ginzu  are  archived  for 
future  reference  and  sent  to  the  user  interface  for  diagnostic  functions  and  user  display.  The 
hardware  used  to  perform  these  functions  is  shown  in  Figure  4-3. 


Figure  4-2:  Local  Ginzu  NILM  Configuration.  Highlighted  blue  area  shows  the  various  functions  performed  following 

raw  data  acquisition.  Figure  courtesy  of  Proper  (5). 
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Figure  4-3:  NILM  Computer  (left)  collects  and  processes  information  received  from  NILM  Sensor  Box  (right)  via 
Ethernet  cable.  Grey  wire  conduit  and  yellow  cable  penetrating  underneath  the  NILM  box  are  raw  voltage  and  current 

measurements  received  from  the  control  panel. 


4.2.1  Detection  Method 

Events  are  detected  when  rapid  power  changes  occur.  Ginzu  uses  a  change  of  mean  filter  to 
perform  this  function.  A  two  column  power  array  of  real  and  reactive  power  is  passed  to  the 
event  detection  algorithm,  where  a  low  pass  filter  is  used  to  remove  the  D.C.  level  of  the 
incoming  data  stream.  The  power  signal  is  then  compared  to  a  user-defined  detection  threshold 
to  determine  whether  or  not  an  event  has  occurred.  Figure  4-4  shows  how  the  detection 
algorithm  receives  the  information  and  produces  event  detection. 


Index  of 
Event 
Detections 


Figure  4-4:  Detector  Change-of-Mean  Filter.  (5) 
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4.2.2  Classification  Techniques 


When  transients  are  detected,  Ginzu  performs  classification  via  two  steps.  First,  steady  state 
power  changes  are  analyzed  in  the  real  envelope  (AP)  and  the  reactive  envelope  (AQ).  Secondly, 
transient  characteristics  defined  by  the  power  spectral  envelope  are  compared  against  a  library  of 
known  system  power  envelopes. 

Steady  State 

When  individual  loads  are  cycled  within  the  system  (i.e.  pump  on/off),  they  produce  a 
corresponding  change  in  the  real  envelope  and  reactive  envelope  and  possibly  other  spectral 
envelopes.  These  changes  can  be  used  as  a  simple  classifier  to  identify  loads.  Ginzu  uses  only 
AP  and  AQ  values  to  classify  events.  Since  the  CHT  system  is  only  comprised  of  four  primary 
loads,  the  need  for  higher  order  harmonics  to  aid  in  classification  was  not  required. 

Transients 

As  described  in  Chapter  Two,  transient-based  recognition  permits  near-real-time  identification  of 
load  operation,  especially  turn-on  events.  Transients  are  identified  by  matching  events  in  the 
incoming  power  stream  to  previously  defined  transient  signatures,  or  “exemplars.”  Pre-training 
has  proven  to  be  a  reasonable  approach  for  very  repeatable  loads  that  show  up  in  large  quantities, 
such  as  fluorescent  lamp  ballasts.  Each  event  detected  is  compared  to  the  full  set  of  exemplars 
by  using  a  least  squares  criterion  to  select  the  appropriate  shifts  and  gains.  The  match  with  the 
lowest  residual  norm  per  number  of  points  is  then  compared  to  a  threshold.  If  the  fit  is  good 
enough,  the  event  is  classified  as  a  match  to  the  exemplar.  If  not,  the  event  is  left  unclassified. 

Correct  classification  of  overlapping  transients  is  possible  using  properly  designed  exemplars  for 
a  system  with  as  few  a  number  of  loads  as  the  CHT  system.  Fingerprint  traces  provide  positive 
identification  of  specific  events  occurring  during  system  operation.  By  comparing  the  shape  of 
the  transients  to  known  system  events,  a  numerical  score  can  be  assigned  to  grade  the  degree  of 
similarity  of  the  two  signals;  this  score  is  known  as  the  correlation  score.  It  is  derived  using  the 
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method  of  least  squares  in  the  Ginzu  algorithm.  This  method  is  discussed  in  detail  in  Lee  (21) 
and  Proper  (5). 

Event  Files 

The  CHT  event  classifier  generates  an  ID  CODE  to  represent  a  specific  system  event  when 
operating  in  its  normal  mode.  TABLE  4-1  lists  normal  system  event  descriptions  along  with 
their  associated  codes.  TABLE  4-2  lists  additional  system  descriptions  that  were  defined  to 
create  a  more  robust  classifier.  Once  a  positive  event  has  been  identified,  an  event  file  is  created 
and  stored  in  the  archive  directory. 

TABLE  4-1:  CHT  CLASSIFICATION  CODES.  (5) 


Description 

ID  CODE 

Vacuum  pump  ON 

(001) 

Discharge  pump  ON 

(003) 

Vacuum  pump  OFF 

(002) 

Discharge  pump  OFF 

(004) 

ALL  pumps  OFF 

(009) 
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TABLE  4-2:  NON-STANDARD  CHT  CLASSIFICATION  CODES.  (5). 


Description 

ID  CODE 

Vacuum  pump  ON,  but  indicates  low  running  power 

(005) 

Vacuum  pump  with  low  running  power  has  turned  OFF 

(006) 

Unclassifiable  event  (probable  noise) 

(000) 

Unidentified  ON  event 

(007) 

Unidentified  OFF  event 

(008) 

Vacuum  pump  cycling  7 

(021) 

Discharge  pump  cycling 

(043) 

Cycle  of  Unidentified  Load 

(087) 

2  Vacuum  pumps  ON 

(Oil) 

Vacuum  pump  ON  and  Discharge  pump  ON 

(013) 

4.2.3  State  Verification 

If  no  transients  are  detected  within  a  given  window,  the  classifier  does  not  need  to  be  called. 
Ginzu  uses  this  opportunity  to  verify  the  current  state  of  operation.  This  is  accomplished  by 
calculating  the  average  power  level  for  a  ten  second  window  and  its  standard  deviation.  These 
values  are  used  in  a  state  verification  function  to  perform  various  checks  and  correct  the  state 
status  if  needed. 


4.2.4  CHT  GinzUI  and  Diagnostic  Package 

The  GinzUI  application  provides  the  interface  between  the  event  file  and  the  NILM  user  (5).  Its 
primary  functions  are: 

•  To  continuously  check  the  user  interface  directory  for  newly  created  event  files. 

•  To  read  event  file  contents  and  move  the  event  files  to  an  archive  directory. 


7  A  ‘Cycling’  load  is  a  load  which  turns  on  but  immediately  turns  off  during  the  same  event  window.  This  occurs 
when  a  template  match  indicates  a  specific  type  of  pump  has  started  -  but  no  relative  power  change  occurs  across 
the  event. 
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•  To  peform  diagnostics  on  event  file  data  and  alert  the  user  if  a  diagnostic  has  failed. 

•  To  allow  the  user  to  graphically  view  event  file  contents. 

Four  CHT  GinzUI  diagnostic  packages  were  available  during  the  September  2008  patrol.  Each 
of  these  was  created  based  on  known  system  faults  that  were  identified  during  previous  NILM 
research  (2)  (20)  A  general  description  of  the  three  faults  is  provided  below8. 

Vacuum  Pump  Clog  Diagnostic 

In  (20),  Piber  describes  a  fault  where  the  vacuum  pump  fails  to  reach  normal  steady  state  power 
level.  This  is  the  result  of  a  clogged  suction  line  or  priming  line  for  the  associated  vacuum 
pump.  This  fault  results  in  abnormally  long  vacuum  pump  run  times  as  well.  Normal  vacuum 
pump  run  duration  is  only  1-2  minutes.  During  this  particular  event,  the  vacuum  pump  remains 
on  for  longer  periods  (several  minutes)  as  it  attempts  to  draw  a  vacuum.  Eventually,  the 
alternating  pump  energizes  to  assist  in  drawing  the  vacuum.  Once  the  high  vacuum  set  point  is 
reached,  both  pumps  are  secured.  Figure  4-5  illustrates  this  sequence  of  events. 


Figure  4-5:  Power  plot  indicating  a  clogged  vacuum  pump  fault.  The  first  pump  energized  runs  for  excessive  periods 
until  the  low  vacuum  set  point  is  reached,  then  the  alternating  pump  assists  until  the  high  vacuum  level  set  point  is 

reached.  (20) 


8  A  detailed  description  of  each  fault  and  its  history  is  provided  in  Piber.  (20) 
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The  clogged  vacuum  pump  is  identified  in  the  classifier  as  ‘VP  ON  with  low  running  power 
(005).  The  diagnostic  function  alarms  when  five  ‘low  power’  vacuum  pump  runs  are  identified 
within  the  ten  most  recent  vacuum  pump  starts.  This  is  accomplished  by  assigning  a  one  to  ‘low 
power’  starts  and  a  zero  to  normal  pump  runs.  If  the  sum  of  the  ten  most  recent  vacuum  pump 
ON  events  reaches  five,  the  CHT  GinzUI  alarms  indicating  a  clogged  vacuum  pump  has  been 
diagnosed.  The  alarm  clears  when  eight  out  of  ten  of  the  most  recent  runs  are  classified  as 
normal. 

Excess  Vacuum  Pump  Runs 

The  vacuum  pumps  are  controlled  via  pressure  switches  located  on  the  side  of  the  sewage 
retention  tank.  One  switch  is  set  to  control  the  operation  of  the  lead  pump  by  energizing  the 
pump  when  pressure  reaches  14”  Hg  and  securing  the  pump  when  pressure  reaches  18”  Hg.  The 
second  switch  controls  the  standby  pump  to  energize  it  if  pressure  reaches  12”  Hg.  If  both 
pumps  are  energized,  the  pressure  switch  de-energizes  both  pumps  upon  reaching  18”  Hg.  (22) 
Each  pressure  switch  has  individual  settings  that  can  be  adjusted  using  dials  within  the  switches. 
If  the  pressure  switch  is  incorrectly  calibrated,  the  vacuum  pumps  may  run  excessively 
attempting  to  reach  a  higher  vacuum  set  point.  This  fault  reveals  similar  power  signatures  as 
those  discussed  for  clogged  vacuum  pumps.  Proper  stated  “The  CHT  UI  continuously  displays 
the  most  recent  vacuum  pump  run  duration  and  the  average  of  the  runtimes  over  the  last  50  runs. 
During  normal  operation,  this  value  should  be  between  30  and  60  seconds.  The  runtimes  are 
calculated  and  retained  as  Vacuum  Pump  OFF  (004)  or  ALL  OFF  (009)  events  (sic).  The  UI 
alerts  the  user  if  the  average  runtime  exceeds  180  seconds.  Since  a  clogged  vacuum  pump 
combined  with  an  improperly  set  or  faulty  sensor  could  result  in  similar  system  behavior,  this 
diagnostic  is  disabled  if  a  Vacuum  Pump  Clog  has  already  been  detected’’. 

Fouled  and  Damaged  Tank  Level  Indicating  Probes 

The  third  diagnostic  algorithm  pertains  to  discharge  pump  operation.  In  March,  2007,  Piber 
observed  erratic  discharge  pump  behavior  based  on  power  plots  of  the  CHT  system.  The 
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discharge  pumps  are  operated  based  on  tank  level  control  circuitry.  This  circuitry  receives  its 
input  via  tank  level  probe  indicators  located  in  the  VCT.  If  these  probes  become  fouled  or  the 
circuit  does  not  function  correctly,  erroneous  discharge  pump  operation  may  occur.  Typical 
discharge  pump  operation  lasts  about  3  to  4  minutes  while  the  pump  removes  120  gallons  of 
sewage  from  the  VCT  to  a  larger  holding  tank.  Figure  4-6  depicts  the  short  burst  associated  with 
a  failed  probe  casualty. 


Figure  4-6:  Failed  probe  casualty  power  response  on  board  the  USCGC  SENECA.  (20) 


The  GinzUI  diagnostic  algorithm  searches  for  this  fault  by  collecting  runtimes  for  each  discharge 
pump  operation.  An  average  healthy  runtime  is  between  3  to  4  minutes.  (5)  If  the  average  of 
the  last  ten  discharge  pumps  is  less  than  10  seconds,  the  GUI  will  alert  the  user  that  there  is  a 
probable  fault  by  changing  the  color  of  the  diagnostic  log  to  yellow  and  indicating  ‘PROBE 
FAILURE’  in  red  (Figure  4-7). 
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Figure  4-7:  USCGC  Escanaba  CHT  GUI  Illustration.  This  screen  capture  indicates  a  possible  probe  failure  based  on 
indicators.  Note  the  diagnostic  log  lists  the  date  and  time  of  the  potential  failure.  (5) 


Leak  Detection 


Vacuum  leaks  in  the  CHT  system  have  been  documented  on  multiple  occasions  during  NILM 
research  (3)  (20)  (5).  Mosman  noted  a  major  vacuum  leak  on  the  USCGC  SENECA  during  his 
studies  (3),  which  led  to  further  studies  on  how  to  determine  if  a  leak  is  present  in  the  sewage 
system.  Mosman  used  statistical  methods  that  modeled  the  time  between  vacuum  pump  runs  as  a 
Poisson  process.  Histograms  were  used  to  visualize  the  distribution  of  the  pump  times.  It  was 
observed  that  major  leaks  in  the  system  were  visible  as  shifts  in  the  histogram  peaks  and  change 
in  the  general  shape.  Piber  followed  this  research  with  field  studies  on  the  ESCANABA  that 
revealed  the  occurrence  of  a  major  leak  using  power  data  captured  by  the  NILM. 


Ginzu  Performs  leak  detection  using  two  separate  methods.  For  larger  leaks,  Ginzu  looks  for 
sharp  edges  in  the  histograms.  If  the  edge  exceeds  a  tuned  threshold,  the  diagnostic  reports  a 
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leak  has  been  detected.  For  smaller  leaks,  Ginzu’ s  leak  diagnostic  looks  for  gradual  changes  of 
the  distribution  over  a  period  of  days.  (5) 


4.3  USCGC  ESCANABA  Local  Ginzu  NILM  Evaluation 

In  September  2008,  the  USCGC  ESCANABA  departed  the  Boston  area  with  two  NILMs 
installed.  A  Ginzu  NILM  was  installed  at  the  local  CHT  control  panel,  while  an  off-line  NILM 
was  placed  further  upstream  to  monitor  aggregate  loads. 9  This  chapter  will  describe  how  the 
local  NILM  provided  the  crew  with  at-sea  diagnostics,  as  well  as  provide  a  critical  review  of 
classifier  and  diagnostic  performance  using  archived  data. 


The  following  quote  was  taken  from  a  letter  of  appreciation  to  the  MIT  LEES  (Laboratory  for 
Electromagnetic  and  Electronic  Systems)  team  on  behalf  of  the  Engineering  Department 
following  the  patrol: 


The  NILM  provided  great  support  to  the  crew  during  our  most  recent  deployment.  With  the 
NILM  mounted  to  the  shipboard  sewage  vacuum  pumps  for  vacuum  collection  tank  (VCT)  we 
were  able  to  monitor  the  run  time  and  frequency  of  the  vacuum  pumps  in  automatic 
mode. .  .Almost  immediately  upon  installation  the  NILM  indicated  that  the  pumps  were  in  fact 
operating  with  a  periodicity  and  duration  that  suggested  excess  operation.  Further  inspection  of 
the  VCT  system  revealed  failing  check  valves  and  pump  deterioration.  Had  the  NILM  not  been 
installed,  the  pump  most  probably  would  have  been  run  to  the  point  of  failure.  The  NILM 
provided  real  time  data  that  allowed  maintainers  to  extrapolate  trend  performance  to  avoid  a 
complete  system  failure  that  would  have  required  emergency  logistic  support.  (23) 


4.3.1  Ship's  Force  Use  of  Diagnostics  to  Identify  Vacuum  Leak 

On  October  19,  2008,  while  deployed,  the  crew  of  the  ESCANABA  replaced  #2  vacuum  pump 
and  motor  assembly  (Figure  4-8)  due  to  excessive  wear  to  the  pump  and  motor  that  was 
identified  during  visual  inspections.  The  cause  of  the  motor  and  pump  degradation  was 
excessive  cycling  of  the  vacuum  pump  due  to  a  major  leak  in  the  system.  It  was  determined  that 


9  The  off-line  NILM  only  gathered  ‘prep’  data.  It  did  not  provide  crew  members  with  a  user  interface  or  diagnostic 
package.  Results  from  the  off-line  NILM  were  discussed  in  Chapter  Two. 
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a  faulty  check  valve  downstream  of  the  vacuum  pump  was  not  seating  properly,  thus  allowing 
high  pressure  air  to  bleed  into  the  vacuum  system.10 


Figure  4-8:  Replacement  pump  for  #2  Vacuum  Pump.  This  photo  was  taken  following  the  ESCANABA’s  return  to  port 

in  November  2008. 


One  indication  of  a  major  leak  in  the  vacuum  system  is  excessive  vacuum  pump  cycling.  During 
routine  system  events  such  as  toilet  flushing,  the  system  pressure  drops  in  a  step-like  fashion. 
Leaks,  depending  upon  their  size,  display  more  gradual  and  steady  decreases  in  system  pressure. 
Mosman  (3)  illustrated  this  in  Figure  4-9  below. 


tm 


Figure  4-9:  Sewage  system  pressure  trace  (upper  plot)  and  time-aligned  vacuum  pump  transients  (lower  plot).  Plots 
illustrate  the  pressure  response  of  the  system  due  to  routine  events  (flushes)  and  slow  leaks  in  the  system. 


10  The  exact  events  leading  up  to  the  replacement  were  not  provided.  The  reconstruction  of  events  is  based  on  post¬ 
patrol  analysis  of  the  data  and  email  correspondence  with  the  ship’s  engineer. 
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As  the  leak  becomes  larger,  the  slope  of  the  system  pressure  increases  and  the  vacuum  pumps’ 
cycling  frequency  increases.  This  system  behavior  prompted  the  crew  to  troubleshoot  and 
identify  the  faulty  check  valve. 


Figure  4-10  shows  the  vacuum  system  behavior  prior  to  and  following  system  repairs.  A  drastic 
drop  in  pump  cycles  per  hour  occurred  after  the  repair,  since  the  system  was  then  capable  of 
maintaining  vacuum  during  no-load  conditions. 


Tirre  (Mnutes) 


Figure  4-10:  Local  power  plots  taken  before  (top)  and  after  (bottom)  repairs.  Each  plot  sample  was  taken  during  the 

0600  hour  for  more  accurate  comparison. 

The  same  information  on  system  cycling  was  provided  to  the  crew  during  the  patrol  via  the  user 
interface  in  a  slightly  different  format.  In  order  to  re-create  what  the  crew  saw  while  deployed, 
one  hour  of  archived  events  (0600-0700)  from  September  20,  2008  was  processed  using  Ginzu. 
The  screen  capture  in  Figure  4-11  displays  the  system  behavior  via  statistics  gathered  by  the  user 
interface  for  the  one  hour  sample.  There  were  26  vacuum  pump  cycles  with  an  average  of  61 
seconds  for  runtime.  The  data  on  the  right-hand  side  of  the  screen  displays  information  on 
vacuum  pumps  and  discharge  pumps.  The  average  discharge  pump  runtime  displays  ‘NaN’, 
because  there  were  no  discharge  pump  events  during  this  one  hour  sample  window.  This 
information  enabled  the  crew  to  note  the  high  vacuum  pump  run  periodicity  and  led  to  corrective 
actions. 
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Figure  4-11:  CHT  GUI  screen  capture  prior  to  repairs. 


Similarly,  one  hour  of  post-repair  data  was  processed  using  Ginzu.  The  screen  capture  in  Figure 
4-12  illustrates  the  results.  The  classifier  correctly  labeled  all  vacuum  pump  events.  The 
average  vacuum  pump  runtime  remained  around  one  minute,  which  is  normal.  However,  the  key 
difference  is  the  number  of  vacuum  pump  runs  per  hour,  which  was  26  prior  to  repairs  and  10 
following  repairs. 


Figure  4-12:  CHT  GUI  display  following  repairs.  Data  was  taken  from  October  10,  2008,  0700. 
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4.3.2  Classifier  and  Diagnostic  Effectiveness 

In  order  to  assess  the  NILM’s  performance  during  the  patrol  period,  archived  data  taken  from  the 
NILM  laptop  was  processed  using  the  ViewerUI.  Ginzu’s  CHT  classifier  was  used  to  analyze 
preprocessed  data  and  generate  event  files.  Once  the  event  files  were  created,  the  ViewerUI  was 
used  to  select  any  event  file  from  a  list  box.  Figure  4-13  provides  an  example  of  a  normal 
vacuum  pump  start  displayed  by  the  ViewerUI. 


Figure  4-13:  A  normal  vacuum  pump  On  event  (001). 


The  viewer  also  has  the  ability  to  simultaneously  display  a  captured  event  and  either  the  vacuum 
pump  or  discharge  pump  template,  which  allows  the  user  to  compare  events  against  their 
exemplars.  The  right-hand  side  lists  the  events  that  are  called  from  the  archive.  Note  that  the 
event  file  highlighted  has  a  classification  code  of  (001),  which  corresponds  to  a  ‘Vacuum  Pump 
ON’  event.  Each  event  window  is  ten  seconds  in  length.  The  green  arrow  along  the  time  axis 
indicates  when  the  event  is  detected;  the  red  arrow  corresponds  to  the  end  of  the  lock-out 
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window.  In  (5),  Proper  explained  that  a  1.25  second  lock-out  window  was  invoked  to  prevent 
multiple  detections  during  a  transient. 

Random  samples  were  chosen  to  analyze  the  overall  effectiveness  of  the  classifier  and  diagnostic 
functions  of  the  Ginzu  NILM.  This  was  only  the  second  time  the  Ginzu  NILM  configuration  was 
installed  for  an  extended  at-sea  period.11  These  random  sample  windows  were  processed  using 
the  CHT  classifier,  and  then  event  files  were  viewed  using  the  YiewerUI.  Eleven  hours  of 
archived  data  taken  from  October  17,  2008  was  analyzed.  A  total  of  168  event  files  were 
captured;  Ginzu  correctly  identified  99%  of  the  events.  The  classifier  was  only  unable  to 
recognize  two  transient  events. 

4.3.3  Vacuum  Pump  Clog  Detection 

The  random  sample  also  revealed  that  the  classifier  function  correctly  identified  clogged  vacuum 
pumps  on  several  occasions.  Recall  from  Chapter  3,  the  classifier  labels  a  potential  clogged 
pump  as  ‘VP  ON  with  low  running  power’.  If  five  out  of  the  last  ten  events  are  classified  as 
such,  the  GinzUI  alarms  to  indicate  a  clogged  vacuum  pump.  The  excessive  vacuum  pump  run 
duration  coupled  with  the  lower  than  normal  steady  state  power  level  is  illustrated  in  Figure  4-14 
using  a  spectral  power  plot.  Figure  4-15  is  the  transient  from  Figure  4-14  that  occurs  nine 
minutes  into  the  hour. 


11  Details  of  first  patrol  period  are  provided  in  Proper’s  thesis  (5). 
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Figure  4-14:  Power  plot  indicating  possible  clogged  vacuum  pump.  The  longer  than  normal  run  time  of  the  first  pump 
coupled  with  the  lower  than  average  steady  power  level  causes  the  Ginzu  classifier  to  classify  event  as  ‘Vacuum  Pump  ON, 

Seal  Clogged’. 


Figure  4-15:  Vacuum  pump  clog  fault.  The  ViewerUI  contains  Ginzu’ s  classification  (top  stamp)  based  on  the  lower  than 

normal  power  and  excessive  run-time. 
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4.3.4  Abnormal  Pump  Transients  and  Multiple  Events 

Data  sets  prior  to  the  pump  and  valve  repairs  were  also  analyzed.  Events  captured  on  September 
29,  during  the  0000-0100  hours  revealed  unusual  pump  signatures  for  one  of  the  vacuum  pumps, 
(b)  Abnormal  pump  start  down  transient  classification 

Figure  4-16  shows  two  consecutive  events  captured  during  this  period.  The  classifier  correctly 
identified  the  vacuum  pump  transient  even  with  the  unusual  up  and  down  transient  that  occurs  as 
the  pump  reaches  its  steady  state  value.  Note  the  trailing  edge  of  the  ‘hump’  is  outside  of  the 
lock-out  window. 
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Figure  4-16:  Classifier  performance  during  abnormal  vacuum  pump  start. 


Figure  4-16  (a)  is  considered  a  ‘multiple  event’  in  a  single  detection  window.  If  multiple  events 
are  determined  based  on  numerical  derivatives  of  the  data  window,  Proper  stated  “...the  event  is 
‘arbitrated’  by  calculating  the  AP  and  AQ  across  the  lockout  window”  (5).  In  this  particular 
case  the  algorithm  functioned  properly  and  classified  the  event  as  a  vacuum  pump  start.  The 
transient  event  (b)  is  the  trailing  edge  of  the  unusual  ‘hump’;  correctly  labeled  as  an 
‘Unidentified  Downpower’  by  Ginzu.  The  power  level  change  due  to  the  trailing  edge  of  the 
hump  did  not  correspond  to  the  values  assigned  for  either  vacuum  pumps  or  discharge  pumps,  so 
it  was  correctly  labeled  as  ‘unidentified’. 

The  next  classified  event  detected  was  when  the  vacuum  pump  was  secured.  This  is  illustrated 
in  Figure  4-17.  Even  with  the  unusual  transient,  the  classifier  eventually  self-corrects  itself 
based  on  power  level  information,  thus  returning  to  an  ‘All  Pumps  Off  state.  Subsequently,  the 
next  event  was  correctly  classified  as  a  vacuum  pump  start  in  Figure  4-18. 
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Figure  4-17:  Classifier  recovery.  The  6 All  Pumps  Off’  event  is  followed  by  the  next  vacuum  pump  transient,  which 

correctly  classified. 


4.4  Performance  Summary 

Overall,  the  performance  of  the  automated  classification  system  was  highly  successful.  The 
Ginz.u  NILM  accurately  detected,  classified,  and  diagnosed  system  transients  throughout  the 
patrol  period  with  no  system  downtime.  With  little  user  training,  the  crew  was  able  to  Ginzu  ’s 
real-time  user  interface  to  diagnose  a  vacuum  leak,  preventing  catastrophic  at-sea  failure. 

The  detection,  classification,  and  diagnostic  methods  were  relatively  simple  schemes  limited  to 
four  operating  pumps.  As  discussed  in  Chapters  Two  and  Three,  as  the  NILM  becomes  more 
‘non-intrusive’,  the  disaggregation  becomes  more  challenging  and  may  require  more  complex 
techniques  for  extracting  target  information.  However,  the  methods  used  for  power  monitoring 
of  the  CHT  system  demonstrate  the  potential  of  the  NILM  as  a  diagnostic  tool. 


93 


Chapter  5 


Conclusions  and  Future  Work 

5.1  Conclusions 

Shipboard  aggregate  power  monitoring  is  a  viable  solution  to  sensor  count  reduction,  which  will 
benefit  the  Navy’s  efforts  in  condition  based  maintenance.  However,  there  are  multiple  issues 
that  challenge  the  extension  of  NILM  in  efforts  to  increase  its  non-intrusive  characteristic. 
Multiple  concerns  were  addressed  throughout  this  thesis.  Three  more  prominent  obstacles  to 
NILM  aggregate  power  monitoring  include:  quantization,  measurement  distortion,  and 
disaggregation  of  loads.  Chapter  2  revealed  that  quantization  and  measurement  distortion  are 
less  likely  to  impact  aggregate  power  monitoring  when  compared  to  disaggregation. 

Effective  load  disaggregation  is  very  application-specific,  thus  various  methods  of  extracting 
loads  to  perform  classification  and  diagnostics  are  necessary  for  successful  aggregate  power 
monitoring.  Specific  methods  of  disaggregation  were  examined  using  ESCANABA  and  LBES 
data.  Results  revealed  that  the  cross-correlation  and  peak  detection  method  was  highly 
successful  on  the  ESCANBA  data.  However,  LBES  data  severely  challenged  this  approach  of 
disaggregation. 

Chapter  4  provided  patrol  results  from  the  ESCANABA  where  the  Ginzu  NILM  was  installed. 
An  application-specific  NILM  designed  to  monitor  the  CHT  system,  the  Ginzu  NILM  performed 
exceptionally  well.  Crew  members  exploited  the  real-time  user  interface  to  diagnose  a  leak  in 
the  CHT  system.  This  real-world  application  of  the  NILM  demonstrates  the  potential  of  NILM 
to  bridge  the  gap  between  data  acquisition  and  real-time  user  information. 
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5.2  Future  Work 

Data  acquisition  improvements  are  currently  in  progress.  The  NerdJack  is  one  promising 
aspect  of  data  acquisition  that  may  lead  to  improved  aggregate  power  monitoring.  Its  increased 
voltage  range  coupled  with  an  increase  number  of  bits  will  yield  higher  data  resolution.  This 
will  enhance  the  NILM’s  ability  to  digitally  “zoom”  into  data  and  look  for  key  features. 
Furthermore,  the  NILM  will  become  more  capable  of  distinguishing  unrelated  transients  from 
target  transients,  increasing  its  ability  to  disaggregate  upstream  data. 

Another  area  of  current  research  that  may  provide  fruitful  results  in  aggregate  data  collection  is 
Zoom  NILM  (24).  The  primary  objective  of  Zoom  NILM  is  to  achieve  18 -bit  resolution  over  a 
wide  range  of  currents.  In  order  to  achieve  this,  various  hardware  changes  are  being  considered. 
For  instance,  the  dynamic  range  of  the  current  transducer  may  be  expanded  using  flux 
cancellation.  Basically,  magnetic  cancellation  of  large  currents  is  used  to  maintain  measurement 
ranges  of  the  current  sensor  within  tolerance. 

As  with  any  research,  field  data  collection  and  analysis  must  continue  in  order  to  test  hypotheses. 
It  was  shown  that  disaggregation  techniques  such  as  cross-correlation  and  peak  detection  work  in 
limited  applications.  Research  throughout  this  thesis  focused  on  one  or  two  attributes  to  perform 
disaggregation.  Future  work  should  also  consider  alternative  solutions  to  extracting  target 
information  from  aggregate  power  streams.  One  possible  solution  is  through  machine  learning 
algorithms  that  group  clusters  of  load  features.  A  cluster  approach  will  provide  additional  details 
of  a  load,  increasing  the  NILM’s  resolution  in  the  aggregate. 


12  A  new  Ethernet-based  ADC  card  being  developed  at  MIT.  Further  description  provided  in  (55). 
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Appendix:  MATLAB  Script  Files 


%SCRIPT  TO  PLOT  DOWNSTREAM  DATA  FROM  ESCANABA  CHT%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%Keith  Douglas %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2' 

clc; clear  all; 

y=textread ( ' snapshot-20081017-060001 . 1 . txt ’ ) ; 

%%%%%%%%%%%%%%%Convert  to  average  power  and  plot%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%RMS  voltage 
V  =  440; 


%Measurement  Resistor  Value  set  with  dip  switches  in  NILM  box 
Rm  =  100; 

%Current  Transducer  Scaling  Factor  for  LEM  LA  55-P 
%  K  =  1/1000; 

%Current  Transducer  Scaling  Factor  for  LEM  LA  205-S  or  LA  150-P 
K  =  1/2000; 

%Voltage  range  across  ADC;  gain  code  =0  corresponds  to  +/-  5v 
g  =  10; 


%ADC  bit  size  or  count  range 
bits=2 A12 ; 


%This  takes  RMS  voltage  along  with  current  input  and  scale  factors 
%and%convert s  to  3  phase  average  power.  The  1/64  accounts  for  prep  FFT  that 
%samples  128  points 

Pl  =  V* (1/64) * (g/bits) * (1/Rm) * (1/K) * (1/sqrt (2) ) *y; 

%Time  vector  converted  into  minutes 

t  =  [0:1: length (Pl)-l] * (1/120) *  (1/60) ; 

%Plot  Figure  in  black  and  white  only 
f igure=  figure ( 'Color' ,  [1  1  1  ] ) ; 

%  Create  plot.  Column  1  corresponds  to  real  power;  Col.  2  Reactive  Power 
plot(t, -Pl(:, 2) ./1000-. 1258, ’Color’, [0  0  0]); 

%  Create  ylabel 

ylabel('Real  Power  (kW) ' , ' FontSize ’ , 12 ) ; 

%  Create  xlabel 

xlabel ( ’ Time  (Minutes ) ' ,  ’ FontSize ’ , 12 ) ; 

%  title ('Local  CHT  NILM'); 
axis  tight; 
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^SCRIPT  to  plot  upstream  data  from  ESCANABA  CHT3 


clc; clear  all; 


%Read  in  the  prep  files  for  both  local  and  upstream  events.  These  files 
%have  already  been  time-aligned  while  in  raw  data  format 
%upstream  power  signal 

x=textread ( ’ snapshot-20081022-020001 . 1 . txt ’ ) ; 


%Convert  to  average  power  and  plot 
V  =  440; 

Rm  =  30; 

K  =  2000; 

g  =  10; 

bits=4096 ; 


%Upstream  Power  value 

%column  1  is  real  power;  column  2  is  reactive 

P=  V* (1/64) * (g/bits) * (1/Rm) *K* (1/sqrt (2) ) *x( : , 1) ; 
t  m  [0:1: length (P) -1] * (1/120) * (1/60) ; 


%Plot  Figure 

f igure=  figure ( 1  Color  1 ,  [ 1  1  1 ]  )  ; 

%upstream  plot 

plot (t, P-868, ’ Color ',[ 0  0  0]); 

title (’Power  Panel  2-2 82-3 ’ , ’ FontSize ’ , 16 ) ; 
ylabel(’Real  Power  (W) ’ ,  ’ FontSize 12 ) ; 

%  ylabel (’ React ive  Power  ( VAR) ’ FontSize ’ , 12 ) ; 
xlabel ( ’ Time  (Minutes ) ’ , ’ FontSize 1 , 12 ) ; 
axis  tight; 
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%SCRIPT  TO  PLOT  DOWNSTREAM  AND  UPSTREAM  DATA  FROM  ESCANABA  CHT%%%%%%%%%%%%%%% 


9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-Q- 


clc; clear  all; 


%Read  in  the  prep  files  for  both  local  and  upstream  events.  These  files 
%have  already  been  time-aligned  while  in  raw  data  format 
%upstream  power  signal 

x=textread ( ' snapshot-20081022-020001 . 1 . txt ’ ) ; 

%downstream 

y=textread ( ' snapshot-20081022-020001 . down . txt ’ ) ; 


%Convert  to 

average  power  and  plot 

V  =  440; 

Rm  =  75; 

%local 

Rml  =  30; 

%upstream 

K  =  2000; 

%local 

K1  =  2000; 

%K  for  downstream  CT 

g  =  10; 

bits=4096 ; 

%downstream  Power  value 

%colsl  and  2  provide  real  and  reactive  coefficients  only 

P=  V* (1/64) * (g/bits) * (1/Rm) *K* ( 1/sqrt (2 ) ) *y ( : , 1) ; 
t  =  [0:1: length (P ) -1 ] * (1/120) *  (1/60) ; 

%upstream  Power  value 

%column  1  is  real  power;  column  2  is  reactive 

Pl  =  V* (1/64) * (g/bits) * (1/Rml) *K1* (1/sqrt (2) ) *x( : , 1) ; 

tl  =  [0:1: length (Pl)-l] * (1/120) *  (1/60)  ; 

%Plot  Figure 

f igure=  figure ( 'Color' ,  [1  1  1 ] )  ; 

%upstream  plot 
subplot (2,1,1) ; 

pi ot ( 1 1 , Pl-86  8 ,  ' Col or’  ,  [ 0  0  0  ]  )  ; 

title ('Power  Panel  2-2 82-3 ',' FontSize ',  16  )  ; 

ylabel('Real  Power  (W) FontSize ', 12 ) ; 

%  ylabel (' React ive  Power  (VAR)  1  ,  1 FontSize 1 , 12 ) ; 
xlabel ( ' Time  (Minutes ) ' , ' FontSize 1 , 12 ) ; 
axis  tight; 

%downstream  plot 
subplot (2,1,2) ; 

plot (t+. 44, -P,  'Color' ,[0  0  0  ]  )  ; 
title ( ' CHT  Control  Panel ',' FontSize ', 16 ) ; 
xlabel ( ' Time  (Minutes )  '  ,  ' FontSize ' , 12 ) ; 
ylabel ('Real  Power  (W) ',' FontSize ', 12 ) ; 
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FontSize 


12)  ; 


%  ylabel (' React ive  Power  (VAR) ’ ,  1 
axis  tight; 
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%Script  for  plotting  3SA  and  L0SP2A,  SLOW  to  OFF  exemplar,  and 
cross%correlation  plots  of  upstream  %and  downstream 
%Keith  Douglas,  April  2009 

%MATLAB  plot  m-file  generation  was  used  to  determine  plot  format  codes 
clear  all;  clc; 

2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2' 
%%%%%%%%%%3SA  DEFINITIONS %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%high  and  low  gain  data  available;  coll  Q,  col2  P 
x=textread ( ’ 3SA_f ullprepsnapshot-2 00806 19-1 10001 . 1 . txt ’ ) ; 

%RMS  VOLTAGE 
V  =  440; 

%bits 
b=4096 ; 

%3SA  has  two  channels-ch . 1  VR=5  ch . 2  VR=20 
%VR  =  20; 

VR=5 ; 

%3SA  has  two  amplification  levels, G,  as  well.  Ch . 1  G— 10;  Ch . 2  G=100; 

G=10 ; 

%G=100 ; 

%KN  is  Fluke  i3000  conversion  ratio  of  .1  mV/A 
k=le-4 ; 

%3  phase  average  power 

P  =  sqrt(3)*V*x(l:431900,2)*VR/(64*b*G*k*sqrt(2) ) ; 

2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2' 

2-2'2'2'2'2'2'2'2'2'2'2'T  ^^7\T  T  ^7X1"^  ■pvTp  Tp  t  1\T  T  T1  T  ^hlVT  C  2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 2- 9- 2- 2- 2- 2- 2- 9- 2- 2- 2- 9- 2- 9- 2- 2- 2- 9- 2- 
ooooooooooooJ_iw0.ttJ_i  J_iwr\L7  i_7llirJ.lM±±±vJl\loooooooooooooooooooooooooooooooooooooooooooo 

%LOSP2A 


y=t extread ( ' LO_f ullprepsnapshot-2 00806 19-1 10002 . 1 . txt ' ) ; 
Rm  =  20; 

k  =  2500;  %K=2500  for  LEM  LA  305P 

VR  =  10; 

PI  =  sqrt ( 3 ) *V*  y ( : , 2 ) *VR*k/ ( 64*b*Rm*sqrt ( 2 ) ) ; 
t  =  [0:1: length (PI) -1] * (1/120) *  (1/60)  ; 


2-9-2-2-2-2-2-2-2-9-2-2-2-2-2-2-2-9-2-2-2-2-2-2-2-9-2-2-2-2-2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2-2-2-9-2-2-2-9-2-2-2-9-2-2-2-9-2-2- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


%%%%%%PLOT  OF  3SA  VS.  VARIOUS  LOCAL  LOADS 


2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'2'9-9-2-9-2-2-9-9-2-9-2-9-9-2-9-9-9-9- 

oooooooooooooooooooooooooooooooooooo 


%  Create  figure 

f igurel  =  figure ( ’ Color ’ , [ 1  1  1 ] ) ; 

%  Create  subplot 
subplotl  = 

subplot (2,1,1, ’ Parent ’ , f igurel, ’ Font Size ’ , 14, ’ Font Name ’ , ’ CambriaMath ’ ) ; 
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%  Uncomment  the  following  line  to  preserve  the  X-limits  of  the  axes 
xlim ( [33  37] ) ; 
box ( 1  on ' ) ; 
hold ( 'all' ) ; 

%  Create  plot  of  3SA 

plot (t-. 004, -P/1000,  ' DisplayName ' ,  ' 3SA ' ,  ' Parent ' , subplot 1 ,  ' Color ',[  0  0  0 ]  ) ; 
%  Create  ylabel 

y label ( { ' Power  (kW) ' } , ' LineWidth ' , 1 , ' Font Size ' , 16 , ' Font Name ' , ' CambriaMath ' ) 
%  Create  legend 

legendl  =  legend ( subplot 1 , ' show ') ; 

set ( legendl ,'YColor',[l  1  1 ] , ' XColor ' , [ 1  1  1 ] , ' Font Size ',14) ; 

%  Create  subplot 
subplot2  = 

subplot (2,1,2, ' Parent ' , f igurel, ' Font Size ' , 14, ' Font Name ' , ' CambriaMath ' ) ; 

%  Uncomment  the  following  line  to  preserve  the  X-limits  of  the  axes 
xlim ( [33  37] ) ; 

box ( ' on ' ) ; 
hold ( 'all' ) ; 

%  Create  plot  of  downstream  load 

plot (t, -Pl/1000, ' DisplayName ',' LOSP  2A ',' Parent ', subplot2 ,' Color ',[ 0  0  0]); 
%  Create  ylabel 

ylabel ( { ' Power  (kW) ' } , ' LineWidth ' , 1 , ' Font Size ' , 16 , ' Font Name ' , ' CambriaMath ' ) 
%  Create  xlabel 

xlabel ( { ' Time  (min ) ' } , ' LineWidth ' , 1 , ' Font Size ' , 16 , ' Font Name ' , ' CambriaMath ' ) 
%  Create  legend 

legend2  =  legend ( subplot2 ,' show ') ; 

set ( legend2 , ' YColor ' , [ 1  1  1 ],' XColor ',[ 1  1  1 ] , ' FontSize ' , 14 ) ; 


-o -o-o -o-o -o-o -o-o -o-o -o-o -o-o -o-o -o-o  Cr  OS  S 


correlation 


ooooooooooooooooooooooooooooooooooooooo 


%Define  exemplar 
% LARGE 

%  exemplar=-Pl (254808:256180)  ; 

% SMALL 

exemplar=-Pl (239050:239100)  ; 

%remove  DC  bias 

exemplar=exemplar-mean (exemplar) ; 

%Create  new  figure 

figure 2  =  figure ( ' Color ' ,  [ 1  1  1 ]  )  ; 

%  Create  axes 

axes 2  =  axes ( ' Parent ' , figure 2 ' , ' FontSize ' , 14, ' Font Name ' , ' CambriaMath ' ) ; 


box ( ’ on ’ ) ; 
hold ( 'all' ) ; 


%  Create  plot  of  exemplar 

t_exemp= [0:1: length ( exemplar ) -1 ] * ( 1 / 12  0 ) ; 

plot (t_exemp, exemplar/1000 ,  ' DisplayName ' , 'OFF  to  SLOW  Exemplar', 'Color', [0  0 
0 ] , ' LineWidth ' , 2 ) ; 
axis  tight; 

%  Create  ylabel 

y label ( { ' Power  (kW) ' } , ' LineWidth ' , 2 , ' Font Size ' , 16 , ' Font Name ' , ' CambriaMath ' ) ; 
%  Create  xlabel 

xlabel ( { ' Time  ( sec ) ' } , ' LineWidth ' , 2 , ' Font Size ' , 16 , ' Font Name ' , ' CambriaMath ' ) ; 


%  Create  legend 

legend2  =  legend ( axes2 ,' show ') ; 

set ( legend2 ,'YColor',[l  1  1 ] ,  ' XColor ' ,  [ 1  1  1],' Font Size '  ,  14 )  ; 


%Cross  correlation  of  exemplar  with  Upstream%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


cross=xcorr (exemplar, -P) / norm (exemplar ) .  A2 ; 
M=length ( cross ) ; 

t_cross= [ -M/2 : M/2-1] * (1/60) * (1/120) ; 


%Create  new  figure 

figure 3  =  figure ( ' Color ' ,  [  1  1  1 ] )  ; 

axes3=axes ( ' Parent ' , f igure3 , ' Font Size ' , 14, ' Font Name ' , ' CambriaMath ' ) ; 
box ( ' on ' ) ; 
hold ( 'all' ) ; 


%  Create  plot  of  cross  correlation 

plot ( -t_cross , cross ,  'DisplayName',  'Upstream  Correlation',  'Color',  [0  0 
0 ] , ' LineWidth ' , 2 ) ; 
axis  tight; 

%  Comment  out  to  view  full  data  file 
xlim ( [33  37] ) ; 

line ([0  60], [.9  .9],' LineS tyle ' , ' — ' , ' LineWidth ' , 1 , 'Color ' , ' k ' ) ; 
line ([0  60], [1.1  1.1],' LineS tyle ' , ' — ' , ' LineWidth ' , 1 , 'Color ' , ' k ' ) ; 

%  Create  ylabel 
ylabel ( { ' Correlation 

Value ' } , ' LineWidth ' , 2 , ' FontSize ' , 16 , ' FontName ' , ' CambriaMath ' ) ; 


%  Create  xlabel 

xlabel ( { ' Time  (min ) ' } , ' LineWidth ' , 2 , ' FontSize ' , 16 , ' FontName ' , ' CambriaMath ' ) ; 


%  %  Create  legend 

legend3  =  legend (' Upstream  Correlation ',' Detection  Window'); 
set ( legend3 ,  ' YColor ' ,  [ 1  1  1 ],' XColor ',[ 1  1  1 ],' FontSize ',  14) ; 

%  Create  plot  of  correlation 
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Color  '  ,  [  0  0 


plot ( -t_cross , cross , ' DisplayName ' , ’Upstream  Correlation' , ' 

0 ] , ' LineWidth ’ ,  2 ) ; 
line ([0  60], [.9  .9],' LineS tyle ’ , ' — ’ , ' LineWidth ’  ,  1 , 'Color ’ , ’ k ' ) ; 
line ([0  60], [1.1  1.1],’ LineS tyle ' , ’ — ’ , ' LineWidth ' , 1 , 'Color ' , '  k  '  )  ; 

%  Create  ylabel 
ylabel ( { ’ Correlation 

Value ’ } , ’ LineWidth ' , 2 , ’ Font Size ' , 16 , ’ Font Name ' , ' CambriaMath ’ ) ; 

%  Create  xlabel 

xlabel ( { ' Time  (min ) ’ } , ’ LineWidth ' , 2 , ' Font Size ' , 16 , ' Font Name ' , ’ CambriaMath ' ) ; 
%  %  %  Create  legend 

legend2  =  legend (' Upstream  Correlation ’,’ Detection  Window'); 
set ( legend2 , ’ YColor ’ , [ 1  1  1 ] , ' XColor ' , [ 1  1  1 ] , ' FontSize ' , 14 ) ; 


%Peak  Detection  Algorithm  Peak  Detection%%%% 


9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'9'S' 


%peakdet.m  file  courtesy  of  Eli  Billauear  at 

%http : / / www . bi 1 lauer .co.il/ peakdet . html  for  %release  to  public 
%domain 


%This  portion  of  the  script  calls  the  peak  detector  function, 

identif ies%where  the  peaks  occur,  then  plots  %the  positive  i.d.’s  (based  on 

threshold%set  by  user) 


[maxtabl , mintabl ] =peakdet ( cross , 0 . 1 , -t_cross ) ; 
k=length (maxtabl ) ; 


%preallocate  matrix 
hold=zeros (k, 2 ) ; 


%look  for  max  values  between  threshold  specified  in  ’if’  statement, 
for  i  =  1 : length (maxtabl ) 

if  maxtabl (i, 2)  <=  1.1  &&  maxtabl (i, 2)  >=  0.9 
hold ( i ) =maxt abl ( i ,  2 )  ; 

else 

end 

end 


%Create  new  figure  and  plot  output  of  maxtab  col  1,  time  from  t_cross,  to 
%see  spikes 

figure 4  =  figure (’ Color ’,[ 1  1  1  ]  )  ; 

axes4=axes ( ’ Parent ’ , f igure4, ’ FontSize ’ , 14, ’ Font Name ’ , ’ CambriaMath ’ ) ; 
box ( ’ on ’ ) ; 
hold ( ’all’ ) ; 

%  Create  plot  of  cross  correlation  results  using  peak  detector  to  find  them 

plot (maxtabl ( : , 1 ) , hold,  ’ DisplayName ’ ,  ’ Upstream  Correlation  Spikes',  'Color ’ ,  [  0 
0  0 ] , ’ LineWidth ' , 2 ) ; 
xlim ( [0  60] ) ; 

line  (  [  0  60],  [.9  .9],  ’Line Style',’  —  ’,' LineWidth ' , 1 ,  ’Color ' ,  ' k ' ) ; 
line ([0  60], [1.1  1.1],' LineS tyle ’ , ’ — ' , ’ LineWidth ’ , 1 , 'Color ’ , ’ k ’ ) ; 
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%  Create  ylabel 
ylabel ( { ' Correlation 

Spikes’ } , ’ LineWidth ',1, ' Font Size ' , 16 ,  ' Font Name ’ , ' Cambr iaMath ' ) ; 

%  Create  xlabel 

xlabel ( { ’ Time  (min ) ’ } ,  ' LineWidth ' , 1 ,  ' Font Size ' , 16 ,  ' Font Name ' , ' Cambr iaMath ' ) ; 
%  %  Create  legend 

legend4  =  legend (  'Upstream  Positive  Detection  Spikes’); 
set ( legend4 ,  ' YColor ' ,  [ 1  1  1 ],’ XColor ’,[ 1  1  1 ],’ FontSize ’,14); 


%Cross  correlation  of  exemplar  with  Downstream 
crossdown=xcorr (exemplar , -PI ) / norm (exemplar ) . A 


9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-9-Q- 


2; 


t_crossdown= [-M/2 : (M-l) /2] * (1/60) * (1/120) ; 


%Create  new  figure 

figure5  =  figure (’ Color ’,[ 1  1  1]); 

axes5=axes ( ’ Parent ’ , f igure5 ,  ’ FontSize ’ , 14, ’ FontName ’ , ’ CambriaMath ’ ) ; 
box (’on’) ; 
hold ( ’all’ ) ; 


%  Create  plot  of  cross  correlation 

plot (-t_crossdown, crossdown,  ’ DisplayName ’ ,  ’Downstream  Correlation’ ,  ’Color’,  [0 
0  0 ] ,  ’ LineWidth  \  2) ; 


%Plot  threshold  window 

line ([0  60], [.9  .9],’ LineS tyle ’ , ’ — ’ , ’ LineWidth ’ , 1 , ’Color ’ , ’ k ’ ) ; 

line (  [  0  60],  [1.1  1.1],  'Line Style’,’  —  ',’ LineWidth ’ , 1 ,  'Color ’ ,  ' k ' ) ; axis  tight ; 

xlim ( [33  37] ) ; 

ylim ( [-1  2] ) ; 

%  Create  ylabel 
ylabel ( { ' Correlation 

Value ' } , ' LineWidth ’ , 1 , ' FontSize ’ , 16 , ’ FontName ’ , ’ CambriaMath ' ) ; 


%  Create  xlabel 

xlabel ( { ’ Time  (min ) ' } , ' LineWidth ’ , 1 , ’ FontSize ’ , 16 , ’ FontName ’ , ’ CambriaMath ’ ) ; 


%  %  Create  legend 

legend5  =  legend (’ Downstream  Correlation ',’ Detection  Window’); 
set ( legend5 ,' YColor ',[ 1  1  1 ],’ XColor ’,[ 1  1  1 ],’ FontSize ’, 14 ) ; 

%Peak  Detector  Algorithm 

[maxtab2 , mintab2 ] =peakdet ( crossdown, 0.1, -t_crossdown ) ; 
k=length (maxtab2 ) ; 

%preallocate  matrix 
ho ld2= zeros (k, 2 ) ; 


for  i  =  1 : length (maxtab2 ) 

if  maxtab2(i,2)  <=  1.1  &&  maxtab2(i,2)  >=  0.9 


109 


hold2 ( i ) =maxt ab2 (  i ,  2  ) ; 

else 

end 

end 

%Create  new  figure 

figure6  =  f igure (' Color ',[ 1  1  1]); 

axes6=axes ( ' Parent ' ,  f igure 6 ,  ' Font Size ’ , 14 , ' Font Name ’ , ' CambriaMath ’ ) 
box ( ' on ' ) ; 
hold ( 'all' ) ; 

%  Create  plot  of  cross  correlation 

plot (maxtab2 ( : , 1 ) , hold2 ,  ' DisplayName ' , ' Downs treamCor relation 
Spikes ',' Color ',[ 0  0  0 ],' LineWidth ', 2 ) ; 
xlim ( [0  60] ) ; 

%  Create  ylabel 
ylabel ( { ' Correlation 

Value ' } r  ' LineWidth ' , 1 , ' FontSize ' , 16 ,  ' FontName ' , ' CambriaMath ' ) ; 

%  Create  xlabel 
xlabel ( { ' Time  (  min 

) ' } ,  ' LineWidth ' , 1 , ' FontSize ' , 16 ,  ' FontName ' , ' CambriaMath ' ) ; 

%  %  Create  legend 

legend6  =  legend (' Downstream  Correlation  Spikes'); 

set ( legend6 ,'YColor',[l  1  1 ] ,  ' XColor ' ,  [ 1  1  1],' FontSize ' , 1 4 ) ; 


